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I . INTRODUCTION 


Research efforts this year have been characterized by 
a strong emphasis in development of new topics while main- 
taining and continuing the refinement and advance of several 
important topics established during the last few years. For 
Jupiter, new modeling capabilities have been developed for 
the peculiar directional features of the lo sodium cloud 
discovered by Pilcher (1980a) and for the extended atomic 
oxygen atmosphere of lo discovered by Brown (1980) . For 
Saturn, preliminary modeling efforts have been initiated 
for the extended ring atmosphere detected by the UV instru- 
ment of the Voyager spacecraft (Broadfoot, et al. 1981) and 
for possible extended atmospheres of some of the smaller 
satellites located in the E-ring. Continuing research for 
Jupiter included a further refinement in the analysis of the 
east-west asymmetry of the lo sodium cloud (see the attached 
preprint; Smyth, 1981a) , additional exploratory modeling of 
the shape of the sodium line profile, and improvement in 
the modeling of the interaction of the lo plasma torus with 
the lo sodium cloud. Continuing research for Saturn consisted 
of refinements in the modeling of Titan's hydrogen torus 
(see the attached reprint; Smyth, 1981b) , while continuing 
research for comets has resulted in the assembly and testing 
of the computer code for the extended coma atmosphere. 
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The resoarch reported for the Jupiter and Saturn systems 
is part of our overall program effort that has been designed 
not only to make full use of the new data provided by the 
Voyager and Pioneer 11 spacecrafts in 1979, 1980 and 1981, 
but also to consider the consistency and relationship of 
this new data to the substantial data base already available 
from ground-based, rocket, and Earth-orbiting-satellite 
observations. For the Jovian planets, model interpretation 
of these data may be used to obtain important information 
about the spatial and temporal distribution of extended 
satellite and ring atmospheres, to characterize and investi- 
gate the interaction of the satellites and rings with each 
other and with the magnetospheres of their parent planets , 
and to aid in understanding the concentration, energy 
density and spatial distribution of charged particles in 
the planetary magnetospheres. 

The research effort reported for comets describes our 
progress in developing a general exospheric model to describe 
both the gaseous and dust coma. The development of this 
model was begun in FY 1980 and is now nearing completion. 
Application of the computer model to cometary data is import- 
ant since it is one of the primary exploratory methods used 
to gain information about the composition, chemistry and 
structure of the comet nucleus. 
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IX. EXTENDED ATMOSPHERES OP 10 


Interest in understanding the atmospheres of the outer 
satellites, their extended natures and their planetary inter- 
actions, has grown significantly in the last decade. Major 
motivation and emphasis of the past few years have been fo- 
cused upon understanding the e>;tended atmospheres of lo, the 
innermost Galilean satellite of Jupiter. To date, three 
different extended neutral atmospheres have been observed 
for lo: (1) sodium, discovered in 1972 by Brown (1974) , 

(2) potassium, discovered later by Trafton (1975) , and 

(3) atomic oxygen, discovered very recently by Brown (1980) . 
Other extended atmospheres may also exist for lo. The exis- 
tence of an atomic hydrogen cloud, initially suggested from 
the measured UV emission near the satellite by the Pioneer 
10 spacecraft (Judge and Carlson, 1974; Carlson and Judge, 
1974; Judge et al., 1976), remains uncertain (Mekler and 
Eviatar, 1980) in the light of the more recently discovered 
UV emissions from the lo plasma torus measured by the Voyager 
spacecrafts (Broadfoot et al. , 1979; Sandel et al. , 1979). 

The detection from Voyager data of gaseous SO 2 on lo (Pearl 
et al., 1979) and the detection from Voyager data (Broad- 
foot et al., 1979; Bridge et al., 1979; Warwick et al. , 1979) 
and Earth-based observations (Kupo et al., 1976; Brown, 1978; 
Pilcher and Morgan, 1979; Pilcher, 1980b; Trauger et al., 
1980) of a plasma torus of sulfur and oxygen ions encircling 
Jupiter, suggest in addition that neutral clouds of SO 2 or 
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its chomical fragments, other than atomic oxygen, may also 
be present in the circumplanetary space. 

Large numbers of observations of lo's sodium cloud have 
been obtained during the last sis s'sn years from ground-based 
measurements and tuis substantial data set is continuing to 
grow. Our past modeling program for lo has, until this year, 
concentrated largely upon studying this data set. Several 
of these studies are important and have been continued. 

The scope of these original sodium studies ' 'has , in the last 
year, however, been broadened to include the larger space- 
craft and ground-based data base now available. This broad- 
ened base of modeling serves to both strengthen our past 
research goals and improve our future capacity to answer 
fundamental questions about the escape of gases from lo and 
the energy balance mechanisms operative in the lo plasma 
torus . 

To provide insight and interpretation of these data 
sets for the extended atmospheres of lo, an active modeling 
program has been pursued at AER during the past two and one 
half years. This modeling effort has, however, been in 
development for about seven years, having its origin at 
Harvard University in 1973 (see Fang, Smyth and McELroy, 

1976; Smyth and McElroy, 1977 and 1978), following the dis- 
coveury of the lo sodium cloud by Brown (1974) . A condensed 
summary of this past and our future projected research pro- 
gram for the extended atmospheres of lo is given in Table 1. 
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A discussion of the |,^Togress made in each area of this re- 
search during the last year is presented below. 

2.1 lo Sodium Cloud 

In the past year, modeD.ing of the sodium cloud has con- 
tinued our past program strategy of emphasizing calculations 
that are exploratory in nature. These exploratory calcula- 
tions have addressed four main interrelated areas: (1) mod- 

eling the two-dimensional intensity pattern of the cloud as 
a function of satellite phase angle, (2) modeling the shape 
of the sodium D-line profiles as a function of sa'‘*«^llite 
phase angle, (3) modeling the interaction of the sodium cloud 
with the lo plasma torus, and (4) modeling the peculiar dir- 
ectional features of the sodium cloud recently discovered 
by Pilcher (1980) and also observed by Goldberg et al. (1980) . 
A condensed summary of the past three-year history and the 
future three-year effort of our research in these four areas 
is presented in Table 1. 

(i) Two-Dimensional Intensity Distribution 

The longest and most active sector of our sodium cloud 
modeling has been in the first area indicated in Table 1, 
that of describing the two-dimensional spatial intensity of 
the region B cloud observed in the sky plane (see Figure 1 
for an illustration of regions A, B, Cl and C2 of the sodium 
cloud) . During the last two years, modeling efforts for the 
two-dimensional spatial intensity of the region B cloud 
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Figure 1: Itegions of the Sodiim Cloud. Region A (the satellite disc seen 

from Earth) , B, Cl arid C2 (the region beyond Cl) are depicted 
in relation to Jrrpiter, lo's orbit, and the reference circles 
of radius 10, 20 and 30 Jupiter radii. The dot-plot orbits 
are for sodium atans ejected from lo in the forward orbital 
direction with velocities between 3 and 20 km sec“l<, 
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luive been focusGd upon understanding certain obf^orvod spatiai 
east-west asynunetries that occur when data are compared for 
satellite orbital phase angles separated by 180 degrees. 

These asymmetries in the absolute intensity (Bergstralh et 
al., 1975 and 1977) as well as the distribution of spatial 
intensity (Goldberg et al., 1978; Goldberg, 1979; Goldberg 
et al., 1981) revealed a small periodic modulation of the 
sodium cloud about the earlier deduced steady state cloud 
description of Smyth and McElroy (1978) . A significant dis- 
covery of our modeling in our previous research year (Smyth, 
1979) was that this periodic modulation could be explained 
by the force extended on sodium cloud atoms as they reson- 
antly scatter sunlight. Since this discovery, the mechan- 
ism has been incorporated in our three-dimensional sodium 
cloud model and calculations have been performed and analy- 
sis applied to verify in considerable detail the earlier 
reported results. The results of this last two years of 
modeling effort is presented in the paper (Smyth, 1981a) , 
for which a preprint is included in the appendix, and will 
not be discussed further here. Future modeling efforts, 
to include the spatial and time dependent sodium cloud sink 
introduced by the oscillation of the lo plasma torus about 
the satellite orbit plane, are also discussed in the pre- 
print. This remains as the only major improvement to the 
model before proper inversion of the east-west asymmetry 
data may be undertaken. 
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(ii) Line Profile Studies 


The second area of investigation for the sodium cloud in 
Table 1 is that of studying the shape of the line profiles 
for D1 and D2 wavelength emissions measured by Earth-based 
observers. Studies of the line profile data are actually 
complementary to studies of the two-dimensional spatial 
intensity data discussed above. This is true since infor- 
mation about different velocity components of the initial 
velocity distribution of atoms emitted by lo may be more 
easily obtained from different data sets. The east-west 
absolute intensity asymmetry data are sensitive to the very 
low velocity components (Smyth, 1980a) , while the line pro- 
file data are sensitive to the high velocity components 
(Smyth and McElroy, 1977; Carlson et al . , 1978; Trafton and 
Macy, 1978; Macy and Trafton, 1980) . The two-dimensional 
sodium cloud image data (which includes the east-west 
intensity distribution asymmetry data) are sensitive to 
the intermediate velocity components (Smyth and McElroy, 

1978) . A balanced approach is required to consistently 
unlock the velocity distribution from these data. 

Our line profile modeling during the past two years 
has been divided into two parts. The first part is to pro- 
vide, as a function of satellite phase angle, consistent 
modeling for both the line profile data and the two-dimen- 
sional intensity data (i.e., the region B cloud) within the 
common emission velocity range where overlap of data exists. 
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Tlio £3Gcond part/ wliich has boon our major emphasis this year/ 

is to understand/ as a function of satellite phase angles/ 

the distribution of higher velocities required/ beyond uhe 

lower-velocity components needed in describing the region 

B cloud/ to accurately reproduce the larger Doppler shifts 

observed in the asymmetric line profile skirts. Model 

results reported last year (Smyth/ 1980a ) had already shown 

that the region B cloud (with initial emission velocities 

-1 

less than 3 km sec ) produces sizeable contributions 
(<200 ma) to the asymmetric skirt for phase angles less than 
70 degrees in the east and less than 250 degrees in the west. 
These contributions result directly from viewing the velocity 
dispersion in the forward portion of the region B cloud 
which/ geometrically/ is somewhat along the line of sight 
for these phase angles (see Figure 1) . For larger phase 
angles/ however/ inherent higher velocity components are 
required in the initial velocity distribution to match the 
observed asymmetric line profile skirt. These inherent 
higher velocity components may be an extension of the 
normal low-velocity ^ inner-hemisphere emission mechanism 
associated with the region B cloud/ or may indicate the 
presence of a collisional-sweeping mechanism driven by the 
Jupiter magnetospheric wind moving past lo. Both mechan- 
isms/ illustrated in Figure 2 , are capable of producing 
the general shape of the line profile for a suitably chosen 
but different velocity distribution. The performing of model 
calculations to investigate these two mechanisms has been 
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our primary emphasis during the past year. Future line 
profile modeling will continue to examine line profile data 
to better understand the nature of the inherent higher vel- 
ocity emission mechanism operative at lo. 

(iii) Interaction with the lo Plasma Torus 
The third area of investigation for the sodium cloud in 
Table 1 is that of describing its interaction with the lo 
plasma torus. The presence of an lo related plasma torus 
was discovered previous to the Voyager 1 encounter with 

Jupiter by Kupo et al. (1976) through ground-based measure- 

. + 

ments of excited state emission from S ions. The presence 
of a much more complex plasma torus, approximately centered 
about the centrifugal equatorial plane with the orbital 
radius of lo, was revealed by the Voyager 1 spacecraft ex- 
periments (Broadfoot et al., 1979; Bridge et al., 1979; 
Warwick et al . , 1979) during its encounter with Jupiter. 

The plasma torus extends radially from about 5 Rj to 7.5 Rj 
and provides a strong sink for sodium atoms, emitted by lo 
at 5.9 Rj, through electron impact ionization. The thermal 
electron component dominates the lifetime process for sodium 
in the plasma torus. Electrons in the outer warm region 

-3 

(>5.7 Rj) , with a number density at 5.9 Rj of order 2000 cm 

5 

and a temperature of about 10 ° K (8.6 ev) , can easily ionize 
sodium, with an ionization potential of 5.14 ev, in a time 
of about one hour. Within the cold inner region (>5.7 Rj) , 
however, the lifetime is significantly larger, having a value 
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of about 5 hours at 5.5 Rj and increasingly very rapidly for 
small radial distances from Jupiter. A spatially asymmetric 
lifetime sink is thus produced in the radial direction near 
lo's orbit by the plasma torus. 

Efforts to incorporate the interaction of the lo plasma 
torus in our sodium cloud model were initiated during the 
past year and a half. To date spatial variations of the 
electron number density within the torus have been described 
in two-dimensions, radially along and normal to the centri-* 
fugal equilibrium plane, with angular symmetry assumed about 
Jupiter. For simplicity the centrifugal and satellite planes 
have thus far been assumed i:o be co-planar. The plasma data 
(Eagenal and Sullivan, 1981) for Jupiter's magnetosphere, 
deduced from the Voyager 1 measurements (Bridge et al., 1979), 
are being used to currently update the plasma properties of 
Bagenal, Sullivan and Siscoe (1980) adopted earlier in our 
sodium model., A significant improvement to the sodium model 
will also result by adopting the more accurate description 
for the spatial variation of the electron energy distribu- 
tion deduced by Scudder and Bridge (1981) . Efforts to in- 
corporate the oscillation of the plasma torus about the 
satellite plane will also be initiated in the coming year. 
These refinements in the description of the plasma torus 
are necessary in order to correctly specify the spatial and 
time dependent ionization history of the cloud atoms in the 
model calculations of the density and D-line intensities of 
the dodi\am cloud. Upon completion of these refinements, 
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definite conclusions about the magnitude of the flux of 
sodium atoms emitted by lo and the location of preferred 
areas of enhanced sodium emission from the satellite's 
exosphere can be deduced by comparison of the model calcu- 
lations and the cloud data. 

(iv) Peculiar Directional Features 

The fourth area of investigation for the sodium cloud in 
Table 1 is that of studying the peculiar directional features 
recently detected by Pilcher (see Hartline, 1980; Pilcher, 
l‘i)80a) . These features, somewhat faint compared to the bulk 
of the region B cloud and usually observed outside of lo ' s 
orbit, are transient in nature and extend at times over 2.5 
Jupiter radii from the satellite. Other feather-like exten- 
sions of the B cloud observed by Goldberg et al. (1980) are 
probably observations of the same phenomena. The features 
observed by Pilcher suggest that localized regions of the 
satellite may eject enhanced high-speed streams of sodium 
atoms with velocities of order 9 km sec 

To facilitate study of these directional sodium features 
a cooperative effort has been established with Pilcher this 
past year. Upon examination of raw data (Pilcher, 1980c) , 
our standard region B sodium cloud model was modified by 
adding a localized jet of enhanced high-speed sodium. The 
intensity of the sodium cloud on the sky plane, as would be 
seen by an Earth-observer, was calculated using this simple 
model for various assumptions about the enhanced sodium 
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stream. Some o£ these results (see Smyth, 1900b; an attached 
copy of the presentation is in the appendix, e ;,imilar in 
character to a recurring transient feature observed by Pil- 
cher and strengthen the hypothesis that enhanced streams of 
sodixim may be emitted by lo. The raw data of Pilcher are 
presently being reduced and will, in the near future, be 
available for direct comparison with our model results. 

2 . 2 lo Potassium Cloud 

The potassium cloud of lo was discovered by Trafton 

O O 

(1975) by observing its emission in the 7665 A and 7699 A 
resonance lines. The detection of potassium emissions was 
confirmed by measurements of Mvinch, Trauger and Roesler 
(1976) and Trauger, Roesler and Miinch (1976) . Later obser- 
vations of the potassium cloud by Trafton (1977) suggested 
that it may be similar to the sodium cloud in many respects, 
but that its intensity is less bright and also more modulated 
by the location of the magnetic equator. More recent data 
for the potassium cloud were reported by Porco, Trauger and 
Karlson (1980) . A significant number of observations (199 
spectra) were very recently reported by Trafton (1980) who 
concluded from preliminary analysis that the potassium cloud 
and sodium cloud are very similiar in character and may be 
ejected from the same region of the satellite by the same 
physical mechanism. 

Modeling of the potassium cloud in the past year has 
not been undertaken because of the lack of available 
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measurements. The recent number of observations reported 
by Trafton (1980) , however, provides the necessary incen- 
tive for developing a potassium cloud model and performing 
interpretive model calculations. This will be the focus 
of our potassium cloud research in the coming year. 

2 . 3 New Gas Clouds of lo 

In addition to sodium and potassium, several other 
neutral gas clouds associated with SO 2 and its chemical 
fragments are likely to exist for lo. To date only an 
atomic oxygen cloud has been detected from ground-based 

O 

observation of its 6300 A emission (Brown, 1980) . Explor- 
atory modeling at AER for the oxygen cloud has already been 
undertaken during the past year. A model for the cloud has 
been developed to calculate its brightness on the sky plane 
for oxygen emission lines , both in the visible and ultra- 
violet, that are excited by electron impact. Preliminary 

O 

modeling results for the visible 6300 A emission of the oxy- 
gen cloud have already been obtained (see Smyth, 1980b; an 
attached copy of the presentation is in the appendix) and 

suggest an atomic oxygen source rate from lo of about 
27 -1 

1-2 X 10 atoms sec to explain the observatxons of 
Brown (1980) . This value of the source rate is, however, 
only preliminary, and may be larger when proper account of 
an appropriate velocity dispersion for the oxygen atoms 
expected from lo is included in the model calculation. 
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The continuation of this exploratory modeling for oxygen 
and its expansion to include atomic sulfur will be supported 
in the future by NASA through the Jupiter Data Analysis 
Program. Exploratory modeling for other, yet undiscovered 
gas clouds of lo, will, however, be incorporated as part of 
our continuing research program here. 



111 . EXl'KWDED ATMOSPIUiUES OF THE SATUl^ SYSl’EM 


Interest in extended atmospheres associated with plane- 
tary rings and satellites began with studies of the Saturn 
system over ten years ago. Early ground-based measurements 
initiated by Franklin and Cook (1969) to search for an atmos- 
phere enveloping Saturn's rings yielded negative results. 

More recent rocket measurements by Weiser, Vitz and Moos 
(1977) have detected a hydrogen Lyman-a emission in the 
vicinity of the planet which may be associated with a hydro- 
gen atmosphere associated with the Saturnian ring system. 

Such an atmosphere could be produced by photo-sputtering of 
watei: ice from the ring material (Carlson, 1980) . Very 
recent measurements of the Saturn system by the UV instru- 
ment of the Voyager 1 spacecraft (Broadfoot et al. 1981) 
confirm the earlier rocket measurements of Weiser, Vitz and 
Moos (1977) and establish the existence of a hydrogen ring 
atmosphere . 

In addition to Saturn's rings, Titan, the largest 
satellite of this planet, has been thought capable of 
generating an extended atmosphere of hydrogen, since it has 
a local atmosphere containing methane (Kuiper, 1944) . Photo- 
dissociation of this atmospheric methane followed by atmos- 
pheric chemistry and finally thermal escape from Titan of 
both atomic and molecular hydrogen (Hunten, 1973, 1974; 

Sagen, 1973; Pollack, 1973) would provide the source, and 
the gravitational capture of this source by Saturn would 
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produce ix Buba tan Uial hydroyea torus about tlio planet 
(Dennefeld, 1973 and 1974; McDonough and Brice, 1973a and 
1973b) . A marginal detection of the Titan hydrogen torus 
was made in 1976 and 1977 and reported by Barker (1977) and 
Barker et al. (1980) from Lyman-a data obtained from the 
Earth-orbiting Copernicus satellite. A positive detection 
of the Titan torus was reported more recently by Judge, Wu, 
and Carlson (1980) from Lyman-a data obtained by the ultra- 
violet photometer aboard the Pioneer 11 spacecraft. Very 
recent confirmation of these two previous measurements has 
now been established from the very recently obtained UV data 
of the Voyager 1 spacecraft (Broadfoot et al., 1981). 

Additional interest in extended atmospheres of the 
Saturn system has been generated by the presence of an en- 
hanced plasma density in the vicinity of the Saturn E-ring 
satellites (Enceladus, Tethys, Dione, and Rhez), which was 
detected by the plasma and uv instruments of Pioneer 11 
(Wolfe et al. 1980; Wu, Judge and Carlson, 1980). The source 
of these ions could result from sputtering or photodissocia- 
tion of ices from these satellites (Wolfe et al. 1980) . In 
this case, a situation similar to lo in the Jupiter system 
may occur for one or more of these E-ring satellites. 

Our past modeliaig program for extended atmospheres in 
the Saturn system, has until the last year, been limited to 
describing the extended hydrogen atmosphere of Titan. Our 
best modeling results, previous to the Voyager 1 encounter 
with Saturn, were recently reported by Smyth (1980b; 1981b, 
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SCO Uho attached reprint in the appendix) . Daring the past 
year, development of a model for Saturn's hydrogen ring 
atmosphere has been initiated and preliminary exploratory 
modeling for possible extended atmosphere of the E-ring 
satellites has begun. The efforts reported here represent 
logical developments of these three areas of investigation. 
Discussion of progress in each of these areas is presented 
below and a condensed summary of how this related to our 
past and future projected research program is included in 
Table 2 . 

3.1 Titan ' s Hydrogen Torus 

The three phases of our past modeling of the Titan 

hydrogen torus are indicated in Table 2. The first phase 
was the development of the model in 1978. This model repre- 
sented an improvement over the earlier model of Fang, Smyth 
and McElroy (1976) by removing their restrictions of a 
circular symmetric torus about Saturn (the limit of long- 
lived collisionless gas clouds) and of a massless satellite 
point source. Results for this new model were first reported 
by Smyth (1978) and were later refined in 1979. With the 
encounter of Pioneer 11 with Saturn in the latter part of 
1979, new data were available to improve the model calcula- 
tions. New model results reflecting the Pioneer 11 encounter 
data were calculated in 1980 and were reported by Smyth 
(1980b; 1981b, see the attached reprint). These model 
results suggest that the escape flux of hydrogen atoms from 
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Observations 



TiUan i» abouU aUomu cin"' »ou ami UuiL the Uorua 

is QuaojitiaXIy cylindrieally symmobric about Saturn, except 
near the satellite where an enhancement in gas density is 
indicated. Future modeling efforts will provide interpreta- 
tive analysis of the Lyman-cx observations recently obtained 
for the 'fitan torus from tho UV instrument of the Voyager 1 
spacecraft (Broadfoot et al., 1981) and similar data antici- 
pated from the encounter of Voyager 2 with Saturn. 


3 . 2 Saturn Ring Atmosphere 

Saturn's rings have been suspected for many years of 
having an atmosphere produced by loss of atoms and molecules 
from surface ices and rocky material. To search for such 
an atmosphere, Franklin and Cook (1969) observed the rings 
when the Barth was near the plane of Saturn's ring in 1966. 
Sunlight reflected from a satellite and passing just over the 
ring surface, was analyzed for absorption or enhancements by 
the ring atmosphere. These measurements yielded negative 
results, providing only an upper limit for the density of 
sodium vapor surx'ounding the ring. Estimates for the bright- 
ness and abundances of hydrogen atoms and OH radicals in a 
ring atmosphere wore made later by Dennefeld (1973, 1974) 
assuming various measurements for release and dissociation 
of surface water ice. Lyman-w brightness for hydrogen of 

O 

1.3 Rayleigh, and 3090 A emission brightness for OH radicals 
of 6.5 Rayleighs were suggested and based upon production 
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ratos of 7.5 x 1023 atoms soc“l and 3x10^^ radicals soc"^, 
respectively. 

Detection from rocket measurements of a 200 Rayleigh 
Lyman-a emission in the vicinity of Saturn was reported by 

Weiser, Vita and Moos (1977) and provided the first evidence 
for a hydrogen ring atmosphere. Carlson (1980) indicated 
that the source rate of about 3x10^^ atoms sec"^, necessary 
to explain this observation, could be produced by the 
mechanism of photo-sputtering of ring water-ice and suggested 
that the dominant loss process for hydrogen was collision 
with the rings. He also estimated the mean number density 
of the ring atmosphere to be about 400 H atoms cm“^, assuming 
the Lyman-a emission was produced by resonance scattering 
under optically thin conditions. Very recent measurements 
of the Saturn system by the UV instrument of the Voyager 1 
spacecraft (Broadfoot et al. 1981) al »o detected the Lyman- 
a emission of the hydrogen ring atmosphere, with brightness 
measurements comparing favorably to those of Weiser, Vita 
and Moos (19 77) . 

To better understand the spatial variation and structure 
of the hydrogen ring atmosphere and the supply of ring gases 
to this atmosphere, preliminary modeling efforts have been 
initiated during the last six months. A model for the 
extended ring atmosphere has been under construction in 
which hydrogen is assumed to be emitted from the ring sur- 
face, with a prescribed radial flux distribution having 
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angular syimnetry around Saturn. This modol is an adapta- 
tion of the Titan torus model discussed above. The initial 
velocity distribution of II atoms at the rings may be pre- 
scribed in the model and will depend upon the ejection 
mechanism operative. Different ejection mechanisms are 
under evaluation. The lifetime of the ejected hydrogen 
atoms is also under evaluation , It may depend upon absorp- 
tion by multi-ring encounters or upon a chemical or ioni'^a- 
tion loss process. Future efforts will be directed toward 
completion of the model development and application of the 
model to ring-atmosphere data. 

3 . 3 E-Ring Satellite Atmospheres 

The possibility of the existence of extended atmospheres 
for the E-ring satellites Tethys and Dione has been raised 
by analysis of plasma data (Wolf et al. 1900; Frank et al. 
1980) and extreme ultraviolet data (Wu, Judge, and Carlson 
1980) obtained in the passage of the Pioneer 11 spacecraft 
through the magnetosphere of Saturn. The plasma data indi- 
cated the presence of an enhanced ion density extending from 
4Rg to 9Rg and approximately bounded by the orbits of 
Enceladus (3.97Rg) and Rhea (8.78Rg). Within this plasma 
torus, two distinct maxima of the ion densities occur, one at 
the orbit of Tethys (4.92Rg) and one at the orbit of Dione 
(6.28Rg). The dominant ion species in the radial intervals 
from the planet of 4. 1-7.1 Rg and 7. 1-8. 2 Rg have been 
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and 0 , i:ani)ocl,LvulY (l‘’i’ank tiL al., .LDUO) . 

Extreme ultraviolet emissions in the radial interval between 
5 and 7 Rg have, in addition, been measured by the UV photo- 
meter of Pioneer 11 and associated with excitation of these 
oxygen ions by electron impact (Wu, Judge and Carlson, 19B0) . 
Beyond the plasma torus, in the radial interval 8,2-16 Rg, 
the dominant ion is II'**. 

The source of atoms for the oxygen plasma toi*us and the 
proton plasma beyond the torus is not known. Since both 
Tethys and Dione have water ice on their surfaces (Gruik- 
shank 1979) and the oxygen ion densities are observed to bo 
maximum at the orbital radius of those two satellites, a 
satellite source is a likely candidate. This suggestion is 
reinforced by the similarity of the overall nature of the 
density and temperature of the oxygon plasma torus in the 
Saturn magnetosphere and the lo plasma torus in the Jupiter 
magnetosphere. Particularly striking in this regard is the 
steep gradient in the ion temperature at the inner portion 
of each torus (Frank et al. 1980; Bagenal and Sullivan 1981) . 
For a satellite source, an extended atmosphere of oxygen or 
OH molecules, provided by dissociation of surface water and 
gravitational escape of its chemical products, could upon 
ionisation (or dissociative ionisation) be the intermediate 
link between the satellite and the magnetosphere. 

Analysis of the extreme ultraviolet data of Pioneer 11 
by Wu, Judge and Carlson (1980) indicated that an ion source 
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rate of about 8 x 10 25 sec"”! would be required to sustain 
the oxygen plasma torus emission rate. If the primary source 
of this oxygen is the satellite Dione, a surface flux of 
0 atoms or OH molecules of about 3x10^ cm"’2sec’"l would be 
required. Such a high flux value is not readily supplied 

g 

by estimated photo-sputtering surface rates of about 1x10 
cm“^sec"l (Carlson 1980 ) , but may possibly result by electron 
impact processes (Wu, Judge, and Carlson 19 80) , although 
laboratory data is required to evaluate this physical process. 
In contrast to the UV data, analysis of the plasma data of 
Pioneer 11 by Prank et al. (1980) suggested a much lower 
plasma torus source rate of 1.5x10^^ to l.SxlO^^ oxygen ions 
sec"l. The analysis of Prank et al. (1980) also favors the 
Saturn rings as the primary source of the torus plasma and 
the satellites, Tethys and Dione, as the secondary sources. 
Preliminary analysis of the newly acquired plasma data from 
the Voyager 1 spacecraft encounter with Saturn (Bridge et 
al., 1981) indicates a magnetospheric structure similar to 
that deduced by Frank et al. (1980) , confirming an enhanced 
density of heavy ions (with a mass-to-charge ratio of 14 to 
16 atomic mass units) in the radial interval from 5 to 7 Rg 
from the planet. Further analysis of Voyager 1 data and 
future analysis of anticipated Voyager 2 data for the 
Saturn system should be very useful in better understanding 
the source of plasma in this planetary magnetosphere. 
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In the past year, efforts to understand the source of 
heavy ions in the plasma torus have been restricted to pre- 
liminary analysis relating the spatial envelopes of extended 
atmospheres for the E-ring satellites to the energetics of 
possible ejection processes. Future efforts for the E-ring 
satellites will focus upon further analysis to evaluate 
(1) the existence and structure of extended satellite atmo- 
sphere, (2) the importance of these extended atmospheres 
as a source of magnetospheric ions, and (3) the initial 
distribution and energy of these satellite ions to the 
planetary magnetosphere. 
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IV. EXTENDED ATMOSPHERES OP COMETS 

Coinatary research represents a broadening and natural 
extension of our previous modeling of extended satellite 
atmospheres and was recently undertaken for NASA in FY 1980 
Large outer satellites and comets actually represent two 
extreme cases of solar system bodies which may have extended 
atmospheres. The large satellites are sufficiently massive 
that the scale height of a light atmospheric gas such as 
hydrogen is nearly equal to the satellite radius, so that 
gravitational escape and formation of an extended atmosphere 
is just thermally possible. Comets, on the other hand, with 
their small size and mass, have enormously large hydrogen 
scale heights compared to their radius, and cannot retain 
molecules released from their nucleus by solar heating when 
near the sun. Comets, therefore, create an atmosphere which 
is by nature always extended. Solar system objects with 
sizes immediate to the large satellites and comets, such as 
Tethys and Dione in the Saturn system, may also have extended 
atmospheres if a source of gaseous material is available from 
their surfaces. 

The extended gaseous atmospheres of comets may be 
divided into two parts: the inner coma and the outer coma. 

Each part of the atmosphere has received considerable theo- 
retical investigation. Hydrodynamic models have been con- 
structed to describe the outward flow in the inner coma, 
since in this region intermolecular collisions dominate. 
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Exosphoric itiodols, based on collision-froc trajectorios of 
gas atoms and molecules, have been developed to describe the 
flow in the outer coma. The thrust of our research effort 
has been to develop improved exospheric models with suffi- 
cient flexibility and accuracy for useful analysis and 
interpretation of cometary data. 

The development of our general exospheric model for 
comets during the last year has involved the testing and care- 
ful assembly of several individually constructed subroutines 
to form the overall desired computer model. This effort 
is part of our overall program objective as summarized in 
Table 2. The model at present is about 90% complete., with 
the remaining 10% of the program development dependent pri- 
marily upon the particular comet atmosphere and the emission 
scenario to be explored. The general atmospheric model is 
based upon exact orbit calculations in three dimensions and 
is suitable for interpretative analysis of both the dust 
and neutral gaseous components of extended cometary atmo- 
spheres. Preliminary test calculations of the assembled 
model are currently in progress for an extended gaseous 
atmosphere of comet Halley. The completed exospheric model 
has been developed to predict the density, column density 
and intensity of the coma for arbitrary viewing geometry 
and as a function of heliocentric distance for given source 
and lifetime conditions. The model should be useful both 
in refinement of analysis of past cometary data and in 
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investigation of new data anticipated in future studies of 
comets. Future efforts will continue implementation of 
model improvements and subsequent application to comet ob- 
servations . 
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ABSTRACT 

A new model Is prcacnicd for Titan’s hydrogen torus, capable of describing its time evolution 
under the influence of the gra^dtalional fields of both the satellite and planet. This description is valid 
when there ore no collisions between cloud atoms. It represents on improvement over earlier models, 
where attention was restricted to a long-lived, collisionlcss gas cloud exhibiting circular symmetry 
about Saturn and to a massless, satellite point source. In tins model, the lifetime of the gas cloud is 
limited by ionization impact of hydrogen atoms witlr the solar wind and planetary mugnetospheric 
plasma and is estimated from Pioneer U data to be of the order of lU^ s. Tor tliis lifetime value, 
measurements of the solar-resonance- scattered La intensity of the torus, obtained both by tire 
Earth-orbitiitg Coperniem satellite and by the UV inslrunKuit of Pioneer li, suggest that 
the hydrogen torus completely surrounds Saturn but that it is not cylindrically symmetric about the 
planet. This lack of synunetry is caused by the finite lifetime of hydrogen atoms in die planetary 
magnctosphcre/solar-wind environment. The time-evolution torus model is thus introduced for data 
interpretation. The measured vertical extent of the torus suggcstii that hydrogen atoms arc emitted 
from Titan’s exosphere with velocities near, and only sUghtiy in excess of, the satellite escape speed, 
Titan’s gravitational field is therefore explicitly included, Modeling results presciUcd provide an 
explanation for the basic features present in both La measurements and suggest a satellite escape 
flux of 1-3X 10** atom.s cm ,s"' and a hydrogen torus containing of the ordijr Of 10’“* atoms. 
Subject headings: plaxiCtS! magnetospheres — planets; satellites — planets; Saturn 


I. INTRODUCTION 

Titan, ll)c largest .satellite of Saturn, is known to have 
an atmosphere containing large detectable concentra- 
tions of CII 4 (Kuiper 1944) and perhaps II, (Trafton 
1972), although more recent measurements (Mdnch, 
Tnaugcr, and Rocslcr 1977) make the latter detection 
less certain. Infrared spectra of Titan’s atmosphere by 
Gillett, Forrest, and Merrill <1973), showing a pro- 
nounced peak at 12 jam, also suggest the presence of 
CjH,,. Photochemical processes in Titan's atmosphere, 
in addition to being able to produce and odicr 
heavy hydrocarbons (Strobcl 1974), may also prosddc 
large escape fluxes of hydrogen atoms and molecules to 
the space around Saturn (Hunten 1973, 1974; Sagan 
1973; Pollack 1973; McDonouglr and Brice 1973), 

McDonough and Brice (1973) first noted that most of 
the hydrogen escaping Titan would be gravitationally 
bound to Saturn and would form a totoidal-shaped gas 
cloud around the planet. They suggested that this hydro- 
gen torus miglu be optically detectable from Earth orbit 
and esdmated its density to be between 1 and 10 ^ atoms 
or molecules cm~\ They also noted that the actual 
density and spatial distribution of the torus were func- 
tions of the Ti Ionian atmospheric escape flux, the radius 
and temperature of the exosphere, and the size and 
strength of Saturn's magnetosphere, which, through 


charged-particle collisions, could provide the dominant 
loss process for hydrogen. McDonough and Brice (1973) 
also noted that satellite recapture of the hydrogen ring 
might provide an additional lo.ss procc.s.s for the torus as 
well as a recycling mechanism to decrease the effective 
net atmosphcric-gas loss rate of Titan, Later work 
(Hunten 1974), however, disputed the recycling idea but 
not the concept of a torus. 

A marginal detection of Titan's hydrogen torus was 
made in 197d and 1977 and reported by Barker (1977) 
and Barker etui. <i9HU) from La data obtained from the 
Earth-orbiting Copernicus satellite. These measurements 
suggest a La brigtuncss of about 200 ±100 raylcii^lis 
near Titan resulting from resonance scattering of sunlight 
by torus hydrogen atoms. Recently, a positive detection 
of the Titan-associated hydrogen cloud was reported by 
Judge, Wu, and Carlson (1980) from Lo data obtained 
by the ultraviolet photometer aboard the Pioneer 11 
spacecraft during its encounter vrith Saturn, The Pioneer 
11 measurements indicated the vertical extent of tlic 
cloud to be about 1.5 Saturn radii above and below 
the satellite plane and indicated an angular extent of the 
cloud along the satellite orbit of at least 5 Saturn radii 
from Titan’s disk. A complete hydrogen torus may exist 
at a weaker signal, but the lixnlted Pioneer 11 data 
preclude a more dcfmitive statement on the angular 
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extent of the cloud. The ob.served angular extent, how- 
ever, indicates that hydrogen is definitely in orbit about 
Saturn since the Lagrange sphere of Titan, the region 
inside which the gravitational field of the satellite 
dominates the planetary gravitational field, has a radius 
of only 0.88 Saturn radii. A complete torus is therefore 
anticipated given the relatively long lifetime estimate of 
10^ s for hydrogen atoms in the planetary environment 
near Titan’s orbit (Judge, Wu, and CarLson 1980), The 
absolute La intensity of the torus measured by Pioneer 
II is difficult to determine because of sensitivity changes 
in the instrument during the encounter, but a pre- 
liminary estimate of 100 rayleighs has been given (Judge, 
Wu, and Carlson 1979). 

In this paper, earlier models for the Titan torus arc 
reviewed in § II. The new model, for Titan’s hydrogen 
torus used to interpret both the limited Copernicus and 
Pioneer II measurements is presented in § III. Improved 
measurements of the hydrogen cloud ore anticipated 
during the encounters of yoyager I and Voyager 2 with 
the Saturn system. The results presented here should be 
helpful in preUminary evaluation of the Voyager data. 

II, EARLIER MODELS FOR THE TITAN TORUS 

Models for Titan’s toroidal gas clouds, based upon 
very simple, collisionlcss, orbital analysis of long-lived 
escaping gases, were first advanced by McDonough and 
Brice (1973) and by Dcnncfcld (1974), These models 
essentially sought to calculate an average toroidal gas 
density /F, 



by estimating the volume K of an angularly symmetric 
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torus about Saturn, the satuliitc emission rate Q (atoms 
or molecules per second), and the lifetime r of the 
orbiting gas in the planetary environment. Results of 
these calculations arc summarized in Table I. The emis- 
sion rate may be alternatively expressed in terms of a 
satellite escape flux (atoms or molecules cm"^ s" '), 

Q=<>S, 

where S is the satellite surface area. 

McDonough and Brice (1973) considered a range of 
emission rates (Trafton 1972) for the escaping gases and 
also different lifetime limiting processes. They estimated 
the toroidal volume to be of the order of the cube of 
Titan’s orbital radius and concluded that gas densities 
were probably in the 1 to 10^ atoms or molecules cm“^ 
range. Dcnnefcld (1974), using an assumed escape speed 
of 1.5 km s " ' to estimate the toroidal volume shown in 
Figure 1 and adopting reasonable lifetime and emission- 
rate values, determined a mean hydrogen-atom density 
of 15 cm"^. For a smaller toroidal volume element 
containing Titan (about 1/15 of that estimated for V), 
Dermefeld (1974) determined the density to be about 3 
times the mean density. Along a line of sight in the 
satellite plane, he then estimated that such an elevated 
hydrogen density would produce an optical depth of 4. 
Assunung, however, that one could see only up to 1 
optical depth, Dcnncfeld (1974) then estimated the La 
brightness to be 95 rayleighs. This conversion of column 
density to optical depth and to solar- resonance- 
scattered La intensity can be seen to correspond to 
hydrogen with an effective gas temperature of 500 K 
from the radiative transfer results of McElroy and Yung 
(1975), which are shown here in Figure 2, properly 
scaled for Saturn’s distance. 


TAULH I 

Models for Titan’s Toroidal Gas Clouds 
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Fki. 1.— Torus Cross sections. Several cross sections of the torus volume in a plane containing Saturn and normal to Titan'.s orbit plane 
: shown, where vertical and radial coordinates are in planetary radii. Denncfcld’s overestimated torus volume 1'' and .smaller adopted 
luntc correspond to the trapezoidaMike cross sections indicated by the heavy dashed lines and solid lines, respectively. For comparison, the 
rrect cro.ss .sections, occupied by atoms emitted isotropically by a moving, massic.s.s. Titan point source, arc calculated by the author and 
dieated by the long-short dashed lines for the two emission speeds of 1.0 and I.S km s " 


An unpublished model for Titan’s torus which treats 
; c-xtreme liigh-dcnsity limit of a long-lived gas cloud 
s developed by Sullivan (1973) using a collision- 
minuted fluid description. Densities in this ease were 
imated between 10^ and 10" cm“^, depending upon 
; a.s.sunied details. The existence of such a dense torus 
unlikely, since heating supplied by absorption of solar 
raviolet below 800 A will cause the torus to decay in 
; order of a year, while the refilling time is of the 
der of 2X10‘‘ years (Huten 1977). Densities around 
^ cm"^ or less are far less objectionable. 



Fig, 2,— Solar-Resonance- Scattering Radiative Transfer 
)dcl for H La Emission. The radiative transfer results for H 
■nputed by McElroy and Yung (1975) arc shown, scaled to 
turn, for an effective gas temperature of 5000 K (.solid curve) 
d 500 K (dashed airve). The optical depth of the gas is indicated 
labeled tick marks on each curve. 


The general problem of long-lived gas clouds emitted 
by satellites was discussed from the orbital mcchaiucs 
perspective by Fang, Smyth, and McElroy (1976). In 
their treatment, the gravitational field of the .satellite 
was ignored, and a model for the cloud density was 
then developed within a two-body orbital mec.hanics 
context. Their paper also presented tlie first, detailed, 
three-dimensional model calculations lor an angularly 
symmetric Titan torus. Normalized, two-dimensional, 
number-density profiles of hydrogen in the cross- 
sectional plane containing Saturn nnd perpendicular to 
the satellite plane were calculated for on effective emis- 
sion temperature of 100 K. Real number densities were 
then obtained by multiplying the normalized number 
densities by the total number of particles in the cloud, 
which requires specification of the lifetime and satellite 
flux of hydrogen atoms. Results for the number density 
and the solar-resonance- scattered. La intensity of the 
torus are given in Table 1 for their adopted values of the 
satellite emission rate and lifetime. The mean density 
values indicated correspond to their extreme outer and 
inner contours. Their La cloud brightness of 1000 
rayleighs is obtained by averaging the corresponding 
column density over a square area, 2 Saturn radii on ft 
side, approximately centered on Titan, and by using the 
radiative transfer results of Figure 2 for an effective gas 
temperature of 5000 K. 

All of the models discussed above have in common 
two simplifying physical assumptions: (1) the toroidal 
cloud is long-lived; and (2) the satellite may be treated 
as a massless point source. The first assumption restricts 
attention to angularly symmetric toroidal gas clouds 
about the central planet. The second assumption simply 
neglects the gravitational field of the satellite and is 
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strictly correct only in the limit of high gas-emission 
temperatures or low satellite moss. For Titan, the second 
assumption is approximately valid if the point source is 
considered to represent the velocity of the emitted gas as 
it o.scapes through the Lagrange sphere of the satellite, 
where the gravitational field of Titan is rougliiy neutral- 
ized by Saturn. 

In order to define more clearly and quantitatively 
what is meant by a long-lived torus, Fang, Smyth, and 
MuHlroy (1976) used the concept of an angular diffusion 
velocity for the gas torus relative to the satellite intro- 
duced earlier by Carlson and Judge (1975). Tliese autl)ors 
noted from simple, two-body, orbital andysis that gases 
escaping from the satellite into a collisionlcss planetary 
environment would diffuse both ahead and bcliind the 
satellite. Gases emitted from die leading face of the 
satellite would have periods larger than the Titan period 
and, consequently, Would eventually trail the satellite 
and populate the space outside of the satellite orbit. 
Gases emitted from the trailing face would have a 
smaller period than Titan and would eventually appear 
to precede Titan and populate the space within the 
satellite's orbit about Saturn. Assuming isotropic emis- 
sion of gas from a satellite point of source with a 
Maxwcll-Boitzmann speed distribution. Fang, Smyth, 
and McElroy (1976) defined an average, forward, angu- 
lar diffusion velocity Sl^jmd an average, backward, 
angular diffusion velocity for the gas cloud. 

For the Titan torus, these averaged, forward and 
backward, angular diffusion velocities are given in Fig- 
ure 3 as a function of the ratio of the effective tempera- 



Fia. 3. — Forward and Backward, Angular Diffusion Velocities 
of Titan's Torus. The average diffusion velocities S2f and arc 
both compared with the simple e.xpression f2, where the ratio of the 
temperature to the atomic weight of the emitted atoms is the 
independent variable. 


turc of the emitted gas to the atomic weight of the 
emitted gas, and arc compared with an approximate 
expression for the average angular diffusion velocity: 


a I rrm / 


Here a is the orbital radius of the satellite, while T is the 
effective temperature, and m is the ma,ss of the atom.s 
emerging from tlic Lagrange sphere with an average 
speed K, relative to the satellite, given by 


\ nm I 


For continuous emission, the forward and backward 
clouds will travel through an angular distance of 360° 
with respect to Titan in a time interval t/Ia of tiie order 
of 


Angular symmetry should be an excellent approxima- 
tion for hydrogen atoms in the Titan torus if die hydro- 
gen lifetime in die Saturn environment is long compMcd 
with either or r^. For an effeedve escape speed K of 
the order of 1 km s~ ' or a temperature 7 of the order of 
200 K for hydrogen atoms, die average, forward and 
backward, angular diffusion times, for the Titan 
toius arc, respeedvely, 3.0X 10° s (or 835 lu) and 2.9X 
10° s (or 807 hr). 

To determine how angularly symmetric the Titan 
hydrogen torus may be, these diffusion times must be 
corKpiued with die photoionization and charge-exchange 
lifedmes of hydrogen atoms operative near the satellite 
orbit. The photoionization lifetime at Saturn’s distance 
is esUmated (Gross and Rasool 1964; McDonougli and 
Brice 1973; Carlson, Matson, and Johnson 1975) to be 
about I-2X 10’ s. TTiis process is, however, dominated 
by charge-exchange of hydrogen atoms with cither solar- 
wind or planetary-magnetospheric protons. 

Analysis prior to the Pioneer 11 encounter with Saturn 
(Siscoe 1978) suggested that Titan’s orbit might be cou- 
tained wholly in the planetary magnetosphere, with an 
estimated corotadng proton density of 1.7 cm"^ provid- 
ing a proton flux of 3.4X 10’ cm“’ s"' and a charge- 
exchange lifedme of 8.4X10° s (2300 hr). Analysis of 
the Pioneer U data (Smith et al. 1980), however, has 
shown that the locadon of the Sunward boundary of the 
magnetosphere is sensidve to the varying, solar-wind 
dynamic pressure and occurred, at times, within the 
orbit of Titan. When the magnetopause is compressed 
within the orbit of Titan, the hydrogen lifetime along 
the unprotected orbit will be characterized by the solar- 
wind plasma having a bulk speed of about 500 km s “ ' 
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and a proton density fluctuating between about 0,2 to 
0.3 cm"'* (Wolfe et al. 1980). This gives a lifetime 
varying between about 5.0-3.3X lO’ s. If, however, the 
orbit of Titan lies wholly within the magnetosphere, the 
hydrogen lifetime is determined by llie corolut>sig plu.smu 
(~200 km s"') with a Pioneer II estimated proton 
density of 0.2— 1.0 cm”’ (Frank 1980) and has a value 
between 7.1-I.4X 10’ s. 

The lifetime of hydrogen atoms in Titan’s torus could 
thus be as short as a few of the diffusion times (r^/j) or 
as long as about 25 diffusion times. For the smaller 
lifetime value (~ 10 ’ s), a complete but angularly un- 
symmetrical torus about .Saturn with on enhanced den- 
sity near the satellite would be produced. Such a torus 
could provide an explanation for both the brighter hy- 
drogen torus observed somewhat parallel to the satellite 
plane by the Copernicus satellite (Barker et al. 1980) 
and the less bright and spatially limited torus observed 
somewhat normal to the satellite orbit by Pioneer 11 
(Judge, Wu, and Carlson 1980). To roughly estimate the 
brightness of sueh a torus viewed parallel to Titan’s 
orbit plane, the above smallest lifetime value of 1.4X 10 ’ 
s is used in the angularly symmetric model results of 
Fang, Smytli, and McElroy (1976) instead of their 
adopted value of 2X lO** s, l^s changes their h;^drogen 
column estimate near Titan from IXIO'** cm"* to 7X 
10 '’ cm”’ and corresponds to a change in the optical 
depth from about 8 to 0.6 using the results of Figure 2 
for an effective gas temperature of 5000 K. Tlicir esti- 
mated La intensity of the Titan torus near the elonga- 
tion point of the satellite orbit is reduced from 1000 
rayleigks to about 190 rayleighs and is thus brought into 
gc^ agreement with the intensity observations of Barker 
et at. (1980). Using the radiative transfer results in 
Figure 2 for the lower effective gas temperature of 500 
K, the La intensity would be reduced to 145 rayleighs. 
The actual or appropriate temperature for the hydrogen 
gas is not known, but depends upon the velocity disper- 
sion that atoms exhibit because of their orbital motion 
in the planetary environment. The effective gas tempera- 
ture produced by this velocity dispersion varies over the 
spatial dimensions of the hydrogen cloud, but is likely 
bracketed by the high and low temperature results of 
Figure 2. 

A lifetime value of 1 X 10’ s corresponds to about 3.4 
diffusion periods for hydrogen aton»s in tlte torus and, 
using the calculations of Fang, Sra'!'}h. and McElroy 
(1976), results in a torus number dattiily in the satellite 
plane ranging from a maximum value of about 125 cm”’ 
at Titan’s orbit to a value of about 2.5 X 10”’ cm"’ at 
radial distances of 5 and 60 planetary radii from Saturn. 
The collisional mean free path of hydrogen atoms in the 
torus, assuming a cross section of 3X 10 ”” cm’ and a 
number density of 125 cm"’, is also about 3.4 Titan 
orbital circumferences. For lifetimes substantially longer 
than 10 ’ s, collisions will therefore begin to be im- 
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portant in the center of the torus. Eventually, the colli- 
sionless assumptions adopted here will no longer be 
valid if the lifetime, for fixed emission conditions, is 
increased sufficiently. The analysis which follows lies 
within the.se permissible bounds. 

in. A NEW MODEL FOR THE TITAN TORUS 

To better understand the spatial distribution of hy- 
drogen in the Titan torus, a model capable of explicitly 
describing its time evolution is required. Such a model is 
introduced here. In addition to removing the long-lived 
time restriction of earlier models, the assumption of a 
massless, satellite point source will also be relaxed by 
explicitly calculating the trajectories of hydrogen atoms 
emitted from the satellite exosphere in the combined 
gravitational fields of Titan and Saturn. The initial 
effects on the motion of hydrogen atoms and also the 
later effects of near encounter or actual atom recapture, 
which arc introduced by Titan’s gravitational field within 
the Lagrange sphere of the satellite, arc therefore ex- 
plicitly included. The model is similar in detail to the 
model developed for the extended gas clouds of lo 
(Smyth and McElroy 1977, 1978) and is based upon 
solving the well-known, circular-restricted, three-body 
equations of motion in three dimensions for emitted 
hydrogen atoms. 

Key input parameters for the model are: (1) the 
lifetime of hydrogen atoms in the Saturn environment; 
(2) the escape flux of hydrogen from the satellite; (3) the 
exosphere radius; and (4) the mean satellite emission 
velocity. A value for the lifetime of the order of 10’ s 
was estimated in § II. A lifetime of 2X 10’ s (5600 hr) is 
adopted here and represents a physically realistic value 
for exploring the amount of angular asymmetry in the 
Titan torus. The escape flux of hydrogen atoms from 
Titan is controlled by chemistry (Strobel 1974) and by 
the processes of verdcal transport operative in the atmo- 
sphere (Hunten 1973). A photochcmical-drivcn conver- 
sion of CH 4 into other heavy hydrocarbons such as 
CiH,,, C 2 H 4 , and C,Hf, will determine the relative 
abundance of H and and suggests (Strobel 1979) an 
escape flux ratio H : in the range 1 ; 9 to 3 : 7. Assum- 

ing a combined escape flux for both H and Hj of about 
10 '® cm"’ s"', based on the satellite surface area, pro- 
vides an upper limit hydrogen-utom escape flux of 3 X 
10 '’ cm”’ s"' and is adopted here for model calcula- 
tions. Estimates of the radius of Titan’s exosphere range 
from 4X10’ to lO"* km (Tabarie 1974), and a nominal 
value of 5000 km is adopted. A mean emission velocity 
of 2.0 km s”' from tUs exosphere is assumed and 
provides the hydrogen atoms with sufficient energy so 
that they escape from the Lagrange sphere of Titan with 
a characteristic velocity of about I km s"'. Hydrogen 
atoms escaping the Lagrange sphere with this velocity 
produce a visible toroidal cloud with a vertical extent 
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Fici. 4,— Tilan H Torus. Column density contours, calculated Fio, 5,— Titan H Torus, Same as Fig. 4 except for a 1 X 10’ s 

from the cloud model for a 2.0 km s~ ' isotropic emission velocity (2800 hr) lifetime, 
and a H-atom lifetime of 2 X lO’ s (5600 hr), are shown as viewed 
from above the satellite plane. 
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Fig. 6,— Titan H Torus. Same as Fig. 4 except for a 5 X 10‘ s 
(1400 hr) lifetime. 


Fig. 7, — Titan H Torus. Same as Fig. 4 except for a 1.‘I4 X 
10® s (400 hr) lifetime. 


comparable to the ±1.5 Saturn radii observed by 
Pioneer II (Judge, Wu, and Carlson 1980). The adopted 
emission velocity of 2 km s ~ ' may therefore be consid- 
ered as an estimated value for the mean velocity associ- 
ated with the yet unknown speed distribution operative 
at Titan’s exosphere. The emission of hydrogen from the 
exosphere is assumed to be uniform and isotropic. More 
complex emission conditions could occur if, for exam- 
ple, there is a strong interaction between the planetary 


magnetosphere and Titan’s atmosphere, but such possi- 
bilities are not considered here. 

For these input parameters, model results showing the 
column density contour plot of the toroidal hydrogen 
cloud of Titan are given in Figure 4 as seen from above 
the satellite plane. The cloud is nearly circularly sym- 
metric, except near the satellite, where the column den- 
sity is peaked. The orbital-angular extent of hydrogen 
from Titan of 5 Saturn radii, observed by Pioneer 11 
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(Judge, Wu, and Carlson 1980), is also indicated in 
Figure 4 and corresponds approximately to tlie en- 
hanced angular column density peak. This suggests that 
the UV photometer aboard Ac Pioneer 11 spacecraft 
may have .seen only the brighter portion of the hydrogen 
cloud near the satellite source. This enhanced column 
density region is more obvious for smaller lifetimes, as 
illustrated in Figure 5, Figure 6, and Figure 7. The 
enhancement will eventually dimini, sh and become part 
of a cylindrically symmetric torus for lifetimes signifi- 
cantly longer than the adopted value of 56(X) hr. 

For tlic 2 km s ' monoenergetiu emission assumed in 
Figure 4 through Figure 7, most atoms escaping the 
Lagrange sphere of the satellite have velocities near 1 km 
s"‘ and produce corresponding, two-body, forward and 
backward, average angular diffusion times of the order 
of 1400 hr. This is larger than the approximate diffusion 
times of 8UU hr obtained earlier, assuming a Maxwell- 
Boltzmann emission velocity distribution. The actual 
average diffusion time will, of course, depend upon the 
exact velocity distribution assumed. These average diffu- 
sion times essentially represent the time it takes for the 
fastest half of the cloud atoms to complete 360° relative 
to Titan. Since the density of the torus is weighted more 
heavily toward the slower diffusing atoms, the apparent 
diffusion time for the column density from Figure 6 is 
seen to be more nearly equal to about 2400 hr. 

The Earth-orbiting Copernicus satellite viewed the 
hydrogen torus nearly parallel to the satellite plane us 
shown in Figure 8, where La intensity contours, assum- 
ing solar-resonance scattering, arc plotted. Model input 
parameters arc identical to those of Figure 4. The 
brightness contours ore also labeled for a reduced hy- 
drogen atom escape flux of IX 10’ cm“^ s“' as indi- 
cated by the values given in parenthesis. These labeled 
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intensities are determined directly from model - 
calculated, column-density values using tlie radiative 
transfer results of Figure 2 for an effective gas tempera- 
ture of 5000 K. The total vertical extent of 3 Saturn 
radii, suggested by Pioneer II UV photometer measure- 
ments, is also indicated in Figure 8. The cloud Intensity 
in the vicinity of the satellite orbit elongation point is in 
the range of 150-400 rayleighs for a flux of 3X10’ 
cm“* s“' and is in good agreement with the upper 
brightness value of 300 rayleighs estimated from the 
Copernicus satellite measurement (Barker ei al, 1980). 
For the reduced flux of 1 X 10’ cm ^ s ', the Intensity 
ranges from about 50-170 rayleighs and is in good 
agreement with the lower Copernicus brightness level of 
100 rayleighs. If the radiative transfer of the hydrogen 
torus is more suitably described in Figure 2 by the lower 
gas temperature of 5(X) K, the La intensity of the cloud 
shown in Figure 8, although somewhat reduced, would 
still be in reasonable agreement with the Copernicus 
observation. In the vicinity of the elongation point of 
the satellite orbit, intensities would be in the range of 
120-240 rayleighs for a flux of 3X 10’ cm“* s“' and in 
the range 45- 135 for the smaller flux of 1X10’ cm“^ 
s"'. The actual effective gas temperature of the Titan 
torus will likely be bracketed by the two temperature 
results of Figure 2. 

During the best Pioneer 11 observation of the near 
Titan cloud, the spacecraft was inside of Titan’s orbit 
and the line of siglit of tlie UV phoiomcter was nonra- 
dial (Judge et al. 1980). The spacecraft was approxi- 
mately 6 Saturn radii from Titan so that the projected 
rectangular field of view of the photometer was about 
0.1 Saturn radii along the satellite orbit and 1.0 Saturn 
radii perpendicular to the orbit. This viewing perspective 
reduces the above column densities and, hence, the 
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TITAN 5600 hr. 



Fig. 8,— Lo Intensity of Titan’s H Torus, Intensity contours, calculated assuming resonance scattering of sunlight, arc shown for the 
satellite plane tilted by 3,25“. An isotropic emission of H atoms with an initial speed of 2.0 km s"', a satellite flux of 3 X 10” cm"^s“* and 
the indicated lifetime were assumed. Brightness contour values in parenthesis arc for a reduced flux of 1 X 10” cm“* s“'. 




-r 


ORIGINAL PAGL 
OF POOR QUALITY 


No. I, 19«| 

TITAN 2800 hr. 


TITAN’S IlYDROGHN TORUS 


351 


•149 

167-1 187-7 / r’^9 



Mox. Signal • 256 Royltighs Max.Signol « 235 Rayleighs 

J•lU, 9, — La Intensity of Titan’s H Torus, See legend to Fig. b. 


observed intensities near the satellite by about a factor 
of 3. Estimated maximum intensities for the Pioneer U 
viewing geometry would then vary between about 130- 
155 rayleighs for a hydrogen flux of 3X lO'^ cm"^ s"' 
and 49-56 rayleighs for a flux of I X lO** cm“^ s“‘, 
using the lower and higher temperature results of Figure 
2v respectively. This is in e,Kcellent agreement witli the 
100 rayleigh estimate given by Judge, Wu, and Carlson 
(1979^ 

In the event that the lifetime of hydrogen atoms near 
Titan's orbit is smaller thari tiie adopted S60U-liour 
value, a more asymmetric or even a partial toroidal 
cloud could exist. Tliis is illustrated in Figure 9, Figure 
10, and Figure 11 for a lifetime of 1 X 10’ s (2800 hr), 
5X10'' s (1400 hr), and 1.44X10'’ s (400 hr), respec- 
tively, for a hydrogen atom flux of 3X lO’’ cm"’ s"'. 
This might occur if the plasma densities are sufficiently 
enhanced by solar wind compression. Labeled contour 
intensities in Figures 9- 1 1 correspond to the higher gas 


temperature radiative transfer results of Figure 2. If. 
however, the atoms live longer than the 5600-hour 
estimate of Figure 8, the column density of the hydrogen 
cloud would increase approximately in proportion to the 
increased lifetime value. The La intensity would in- 
crease less rapidly. Near the elongation point of the 
satellite orbit, where, for example, the optical depth in 
Figure 8 is already between I and 2 for the 3X10“' cm 
s '■ • flux vrdue and an effective gas temperature of 5000 
K, a tenfold increase in the column density would only 
increase the intensity to values slightly in excess of 1000 
rayleighs. The model calculation would at this point, 
however, cease to be valid since collisional redistribution 
of hydrogen in the torus would then be important. 

Eventual loss of hydrogen atoms from the torus 
through ionization may provide a significant ion source 
for the outer magnetosphere. An average radial distribu- 
tion of hydrogen atoms in the torus may be calculated 
by integrating over the angular and vertical dimensions 
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Fig. 10,— La Intensity of Titan’s H Torus. See legend to Fig. 8. 
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Max. Signal ■ 149 Rayleighs 

FiCJ. 1 l.—La Intensity of Titan's H Torus. Sec legend to Fig, 8. 


of the cloud. Dividing this distribution by the lifetime of 
cloud atoms provides an average radial distribution for 
the loss rate of hydrogen atoms (given in Fig. 12). 
Assuming a steady-state cloud this hydrogen-atom loss 
rate is the rate at which new protons arc introduced into 
the rotating planetary magnetosphere. Since the intro- 
duction of a new proton requires the loss of an old 
proton, for the charge-exchange lifetime process domi- 
nant in Titan’s torus, only the energy distribution, not 
the number of protons is changed. Most of the neutral 
hydrogen atoms released from the charge-exchange pro- 
cess are on gravitational escape orbits from the planet, 
A net production of protons will, however, result from 
electron impact ionization of torus hydorgen atoms, 
with a rate of approximately 10"^ times the charge- 


exchange rate, and also from photoionization, Such a 
signature in either the energy or the abundance of 
protons ncfir Titan’s orbit may be present in the Pioneer 
11 data or in anticipated Voyager ion data 

IV. CONCLUDING REMARKS 

Observations of Titan’s hydrogen torus and earlier 
modeling efforts are reviewed. Estimated lifetimes for 
hydrogen atoms near Titan’s orbit of the order of 10’ s 
based upon recent Pioneer 11 measurements suggest that 
the torus completely encircles Saturn and is angularly 
unsymmetric, having an enhanced gas density near the 
satellite. New model calculations presented confirm tliis 
hypothesis and, furthermore, provide an explanation for 



Fig, 12.— Torus H-Atom Lo,ss Rate. Model results for the average radial distribution of the loss rate of H atoms from the Titan torus by 
charge-exehangc processes are shown. The total atom population and overall loss rate are 4.6 X lO’'^ atoms and 2.3 X lO” atoms s”', 
respectively. Model parameters are the same as those adopted for Fig. 8, Radial displacement is in units of Saturn radii. 
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ilw loras Ueiccictl hoih by il»e Earih-orbiiinu Copenmm 
sutuliiii; und die UV instrument of Pioneer II, Agree* 
ment between calculated and observed La intensities 
suggests a hydrogen escape flux between I X 10'* cm"* 
.s ■ ' and 3X 10'' cm '* .s"' .should be operative at Titan. 
Tliis produces a torus containing some 10*^ hydrogen 
utom.s. 

The average rate at which hydrogen atoms are lost 
from the (onis i.s calculated. Thi.s loss occur, s mainly 
tljrougli charge-exchange, which modifies die energy 
distribution but not the abundance of the planetary 
ions. A radial profile for this ion injeedon rate is pre- 
.sented. Loss of hydrogen atoms through electron-impact 
ionizaUon will, however, provide a net source of magne- 
tosphcric protons, as will photoionization. The actual 
spatial distribution of satellite magnetospherie ions will, 
of course, depend upon the balance established between 
the ion mjcclion rale of die Titan torus und die ion loss 
rates provided both by ion diffusion inside the magneto- 
sphere und solar-wind ion capture outside of the magnc- 


lo.sphere. The capture lo.s.s rate may be .sub.siantial .since 
die solar wind/magnetospherc boundary is variable and 
can occur at times within Titan’s orbit. 

Tlie Titan hydrogen torus model presented should be 
useful in preliminary evaluation of both the UV and 
chorgcd-particle data anticipated from the Voyupier mis- 
sion during the next year. I'uture modeling of these new 
measurements should moke it possible to refine our 
knowledge of the local Titan atmo.splicrc a.s well a,s of 
the supply and diffusion of ions in the planetary magne- 
tosphere, 

The author would like to thank R, Goody und D. M. 
Hunten for helpful comments. This research was sup- 
ported by the Planetary Aimo.sphercs program of the 
Nudonul Aeronaudes and Space Aduiiziistradon under 
grant NASW-3I74. Acknowledgement is al.so made to 
the Nadonal Center for Atmospheric Research, wliich is 
sponsored by the National Science Foundation, for the 
computing time used in model computations. 
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Introduction 

Recent advances made in modeling Titan's hydrogen torus , 
lo's sodium cloud and lo's oxygen cloud will be presented. 
Model rf:>jsults for these neutral satellite gas clouds 
are based upon following the trajectories of many atom orbits 
emitted by the satellite. 

This involves solving the equations of motion of each 
atom in three-dimensions, for prescribed initial conditions, 
and including proper weight factors along each orbit to simu- 
late both ionization and excitation phenomena. 

Titan's Hydrogen Torus 

A model calculation for Titan's hydrogen torus is shown 
in the first two slides. 

On the left, Lyman-a intensity contours for the hydrogen 
torus, produced by resonance scatte :ing of sunlight, are 
shown as they would appear to the Earth-orbiting Copernicus 
satellite. 

On the right, column density contours are shown as viewed 
normal to the satellite plane, revealing a complete torus 
about Saturn. 

Values assumed for model parameters are indicated in 
the box. 

The density of the torus is limited by the lifetime of 
H atoms in the Saturn environment, here estimated from Pion- 

7 

eer 11 plasma data to be about 2x10 sec. 

The Lyman-a intensity of 200 ±100 Rayleighs, observed 
near the elongation points of Titan's orbit by the Copernicus 
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satellite, is in agreement with the model calculation if 

. 9 - 9 

the satellite escape flux xs between 1x10 and 3x10 atoms 
cm sec 

The angular extent of the hydrogen torus about Titan of 
+5 Saturn radii, observed by the UV photometer of Pioneer 11, 
is indicated in the right slide and corresponds to the en- 
hanced hydrogen density near the satellite source. 

•» 

lo Oxygen Torus 

A model calculation for the lo atomic oxygen torus, 
recently discovered by Brown, is shown in the next two slides 
The viewing perspective shows lo as it appeared at the 
mid-point of Brown's measurements. 

On the left, brightness contours calculated for the 

O 

6300 A oxygen line, assuming electron impact excitation, are 
indicated in Rayleighs and assume a satellite escape flux of 
10 atoms cm sec . 

On the right, column density contours are shown as 
viewed normal to the satellite plane, revealing at the 
torus occurs mainly ahead of the satellite and inside of 
lo's orbit. 

This shape is determined primarily by the lo plasma 
torus, which provides a spatially non-uniform sink for atomic 
oxygen, through electron impact ionization. 

Plasma conditions assumed here, reflect the Voyager X 
encounter conditions, and were recently calculated by Sheman- 
sky using both Voyager UVS and plasma data. 
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The two locations of Brown's observing aperture are 
indicated by the shaded rectangles, and his average measured 
value of 8 ±4 Rayleighs is in agreement with this model cal- 
culaton, if a slightly reduced satellite escape flux of 

Q -2-1 

5x10 oxyge.;’^ atoms cm sec is adopted. 

For the cooler plasma torus conditions of the Voyager 2 
encounter, the model calculated oxygen densities are larger 

•t 

by about a factor of two. 

Under these conditions, the satellite escape flux, 

necessary for agreement with Brown's observation, is reduced 

9 -2-1 

to about 3x10 atoms cm sec 

Id's Sodium Cloud 

Faint linear-like features of the lo sodium cloud, 
recently detected by Pilcher, suggest that localized regions 
of the satellite may eject enhanced high-speed streams of 
sodium atoms. 

Preliminary modeling to explore this suggestion xs 
presented . 

The next two slides show column density contours of the 
sodium cloud viewed normal to the satellite plane for two 
different model calculations. 

These results were obtained by adding to our standard 
region B sodium cloud model, a localized jet 10® x 10® in 
size and centered 20® above the satellite plane. 

The jet is emitted radially from lo with a speed of 9 
km sec ^ and its vector direction, projected onto the satel- 
lite plane, is indicated on each slide by the heavy arrow. 
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Viawocl from tho Earth, with Xo noar the elongation point 
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of its orbit, tho brighter portions of these two jets look 
somewhat alike and are similar to a recurring transient 
feature observed by Pilcher. 

However, if Xo is 45“ ahead of its western elongation 
point, the two jet features may easily be distinguished by 
an Earth observer a« illustrated in the last two slides. 

Future" comparison of data and model results are required 
to understand and hopefully identify localized sodium jet 
emissions . 
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Sodium Cloud with Jet Feature 
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Sodium Cloud with Jet Feature 
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Abstract 


A new model has been developed for the sodium cloud of 
lo capable of describing its phase-angle dependent interaction 
with the solar radiation pressure arising from resonance 
scattering of sunlight in the and D 2 lines. This model was 
developed to more quantitatively explore the hypothesis, ori- 
ginally presented in Paper I (Smyth, 1979) , that solar ra- 
diation pressure provides a mechanism for explaining the east- 
west cloud distribution asymmetry discovered by Goldberg et 
al. (1978) and the east-west cloud intensity asymmetry dis- • 
covered by Bergstralh et al. (1975) . Model calculations 
presented here confirm the original hypothesis and uncover 
explicit mechanisms responsible for the two observed east- 
west asymmetries. Model results presented for the cloud in- 
tensity did not include the spatially and time dependent 
sodium lifetime produced by o£iQi|,,Xation of the lo plasma torus 
about the orbital plane of the satellite. Modifications to 
these results anticipated upon inclusion of the lo p-Xasma 
torus in the model were, however, discussed and represent 
the only major improvement to the model required before proper 
inversion of the east-west asymmetry data may be undertaken. 
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1. Introduction 


In an -earlier paper (Smyth, 1979; hereinafter referred 
to as Paper I) , an explanation of the observed east-west 
asymmetries of the lo Sodium cloud, based upon the effects of 
solar radiation pressure, was presented. These asymmetries 
were discovered by Bergstr.alh et al, (1975, 1977) and 
Goldberg et al, (1978) when Earth-based observations of the 
D-line emissions of the sodi\am cloud were compared 'for dia- 
metrically opposite satellite phase angles east and west of 
Jupiter. In Paper I, the force experienced by sodium atoms 
as they resonantly scatter sunlight in the D-lines and the 
ability of this force to alter the cloud atom orbits were 
explored. Orbit calculations presented showed that this 
solar radiation pressure has two major effects on the cloud. 
First, it acts in a time-dependent fashion, compressing the 
cloud near eastern elongation and expanding the cloud near 
western elongation, thereby providing the correct behavior 
to explain the east-west intensity asymmetry observed by 
Bergstralh et al. (1975, 1977). Second, it causes the 
sodium cloud to tilt closer to Jupiter when lo is near western 
elongation, thereby providing an explanation for the east- 
west cloud distribtuion asymmetry observed by Goldberg et al. 
(1978) . An illustration from Paper I showing these two effects 
is given in Figure 1. 

The major objective of Paper I was to simply demonstrate, 
using orbit calculations restricted to the satellite plane. 


3 



bhat solar radiation pressure provides a mechanism to quali- 
tatively explain the observed east-west asymmetry features . 

The object of this second paper is to further test this me- 
chanism by providing quantitative and more realistic three- 
dimensional model calculations for the sodium cloud. To 
accomplish this, a sodium cloud model, including the inherent 
time dependent effects of solar radiation pressure on the 
cloud atom orbits, has been developed. This model represents 
a substantial modification of the earlier sodiiim cloud 
model of Smyth and McElroy (1977, 1978) and is presented in 
Section 2. 

In addition to solar radiation pressure, the actual in- 
tensity distribution of the observed sodium cloud will depend 
upon the initial velocity dispersion of sodium ejected from 
lo and also upon the spatial and time variations of the sodium 
lifetime in the Jovian environment. Knowledge of this velo- 
city dispersion should be obtained from a careful inversion of 
the sodium data using the above described cloud model into 
which the lifetime information had been properly incorporated. 
The emphasis of the three dimensional modeling presented in 
this paper will hot, however, be data inversion. It is rather 
to provide the preliminary and necessary step of documenting 
the physical changes that occur in the sodium cloud intensity 
because of solar radiation pressure. In this documentation, 
the effects of solar radiation pressure on those 
components of velocity that would make up the actual initial 
/ velocity dispersion of the ejected sodium atoms and that would 
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also contribute dominantly to the cloud int^::r^ity are inves- 
tigated. Changes in the cloud density, introdw'ed because 
sodium atoms arc lost by ionization in a spatially non-uniform 
manner through electron impact with the lo plasma torus, have 
not been included. These changes represent a further numerical 
tailoring of the cloud density and D-line intensities within 
the three-dimensional sodium envelope which has already been 
altered by the action of solar radiation pressure. It is 
the altered intensity distribtuion within this sodium envel- 
ope, without this additional numerical tailoring, that is 
the subject of this paper. The additional effects of this 
tailoring are currently under evaluation and will be con- 
sidered in a future publication. A brief discussion of the 
anticipated effects of this tailoring is, however, included 
in the following section. 
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2. Sodium Cloud iModel 


An improved model for the lo t-odium cloud is presented 
for calculation of its density and solar resonance scattered 
D-line intensities. The model is based upon following the 
trajectories of many atom orbits emitted by the satellite, 
subject to the combined gravitational fields of Jupiter and 
lo and the acceleration produced by resonance scattering in 
the D-lines. This involves numerically solving the equations 
of motion of each atom in three-dimensions, for prescribed 
initial conditions, and including proper weight factors along 
each orbit to simulate both ionization and excitation pheno- 
mena. A description of the equations of motion and the 
resulting sodium cloud model, when the acceleration produced 
by solar resonance scattering is not included, was previously 
considered in detail by Smyth and McElroy (1977, 1978). 
Similar approaches have also been adopted in other sodium 
cloud modeling efforts (Matson et al., 1978; Goldberg et al. , 
1980). In these earlier papers, the motion of sodium atoms 
in the lo-Jupiter system required numerical solution of the 
circular restricted three-body problem in three diraensions. 
The present account pro^^ides an extension of these earlier 
papers by explicitly including the acceleration of sodium 
atoms produced by solar resonance scattering in the D-lines. 

The trajectories of sodium atoms in the improved model 
are governed by what will be called the modified circular 
restricted three-body equations of motion, given as follows; 
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with 


t\ * A A 

X - 2y = + b • (i^cos t + i^sin t) 
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(9) 

(10) 


where b is the acceleration experienced by sodium atoms 
because of resonance scattering in the D-lines, expressed in 
dimensionless form. The normal circular restricted three-body 
equations of motion are recovered when b is equal to zero 
(see Smyth and McElroy, 1977) . 
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The coordinates (x,y,z) measure the displacement of the 
sodium atoms with respect to the center of mass of the Ju- 
piter-Io system in units of the constant separation distance 
between the planet and satellite. Velocities (x,y,z) , accel- 
erations (x,y/Z) and time t are defined in dimensionless 
fashion using the reciprocal of the angular frequency for 
motion of the planet and satellite around their Center of 
mass as a characteristic time. Satellite and planet masses 
are denoted by mi and mg respectively. 

The coordinates (x,y,z) are defined in a non-inertia 
frame which rotates in such a manner to maintain the planet 
and satellite at fixed positions on the x-axis as shown in 
Figure 2. These rotating coordinates are related to a 
dimensionless inertia frame coordinate iK tZ) t with right 
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handed unit vectors (i- ,i ,i^) , by the transformation 
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The inertia frame, also with origin at the center of mass 
of the Jupiter-Io system, is positioned with its C~axis normal 
to the satellite orbit plane and its ?-axis in the plane 
defined by the Earth-Jupiter line and the ^-axis. The 
direction of positive 5 is away from the Earth toward Jupiter 
when the tilt angle of the satellite plane view from the Earth 

A 

is zero. The units vector i is directed away from the 

sun 

sun along the Sun- Jupiter line. 
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ORIGINAL PAGE IS 
OF POOR QUALITY 


The acceleration b is given by 


(bj + bj, )/ 


G<mi + m ) 


-] 


( 12 ) 


where b and b are the contributions to the acceleration 
1 2 

of a sodium atom by resonance scattering in the and 
lines respectively^ where G is the gravitational constant 
and where a is the constant separation distance between 
the planet and satellite. The acceleration b is made 
nondimensional by the characteristic acceleration given in 
brackets, which is effectively the gravitational acceleration 
experienced by a sodium atom in circular orbit of radius a 
about Jupiter. The magnitude of the acceleration ’ b, as 
discussed in detail in Paper I, is as large as 0,01 to 0.02 
and depends upon the radial velocity v of the sodium atom 
with respect to the sun in the following manner: 


hv . 

i j. 

m^,c 1 


(13) 


where 


Yi(v) 



(14) 


Here ^ is the momentvim change experienced by a sodium 

c 

atom of mass m^ upon absorbing a photon of energy hv^^ 
in the line (i=l) or line (i^2) and is the corres- 
ponding rate at which photons are absorbed. In the ex- 
pression for R is the Sun-Jupiter distance in astronom- 
ical units, f- is the oscillator strength of the i line, 

X ^ 
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V 


ORIGINAL PAGE IS 
OF POOR QUALITY 

ifF^^ is the solar continuum- level photon energy flux between 
the and lines at R = 1 , and (v) is the fraction of 
this solar continuum flux available to the sodium atom as 
it is Doppler shifted out of the bottom of the solar D-line 
Fraimhofer absorption feature by the instantaneous radial 
velocity v of the sodium atom relative to the Sun. 

The velocity v is the sum of the radial component of the 
velocity of the center of mass of the Jupiter-Io system 
u^j^ and the velocity of the sodium atom u^^Qjj^in the inertia 
frame with respect to the Sun: 
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<^cm ^atan 
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sun 


(15) 


where 
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atom 
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and where (|,fi/?) is given in terms of (x,y,z) and (x,y,z) 
by differentiation of (11) with respect to t; 
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The presence of the solar radiation acceleration in 
the equation of motion (1) - (3) makes the differential 
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equations explicitly time dependent in two ways. First there 
is the obvious explicit time dependence in the (x,y.,z) frame of 
the changing vector direction of b indicated in equations (l)-,(3) . 
Second there is the complex time- dependent modulation of the 
magnitude of the acceleration vector b as described by 
(13) - (17) . The value of b varies more than an order of 
magnitude and is strongly dependent upon the satellite phase 
angle through the first term in (17) and is also dependent 
upon the instantaneous velocity of the sodium atom (x,y/Z) rela- 
tive to lo through the second term of (17) . These two sources 
of explicit time dependence severely complicate the solutions 
of the equations of motion (1) - (3) . The natural circular 
symmetry of the gravitational problem, that exists in the 
(x,y,z) frame when b=0, is destroyed, thu.s causing each atom 
orbit, ejected from lo with the same initial conditions, to 
have a different trajectory for different initial locations 
of the satellite on its circular orbit. This breakdown of 
the natural circular symmetry of the gravitational problem 
caused by the introduction of the acceleration b gives 
rise to east-west differences in the cloud atom orbits 
dicussed in Paper I and provides a mechanism to explain the 
east-west asymmetries observed in the sodium cloud. 

Model calculations of the density and intensity of the so- 
dium cloud are based upon determining the locus of cloud atoms 
in space and their instantaneous velocities resulting from con- 
tinuous ejection of sodium from lo. This locus of atom positions 


; -^r - . 
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for each individual initial conditions used to describe this 


ejection process, cannot be obtained by simply integrating 
the equations (1) - (3) once over the effective lifetime of 
the atom in the Jovian environment for one orbital location 
of lo as was accomplished in earlier modeling efforts where 
b is zero. In the present case, as the satellite is moving 
around its orbit continuously emitting sodium, the atoms 
ejected at different satellite orbital locations experience 
a different force and therefore have different trajectories 
even though their initial conditions relative to the satellite 
are identical. The determination of the desired locus of 
atom positions in the (x,y,z,) coordinates and their corres- 
ponding velocities therefore must be constructed, for each 
initial condition, from the solutions of the equations of 
motion (1) - (3) for a succession of different locations of 
the satellite along its orbit. The length of the angular 
section, over which these successive solutions are required, 
is determined by the effective lifetime of the sodium atoms, 
and this angular section must be properly pre-positioned 
along the satellite orbit to correctly define the cloud his- 
tory at the time of observation. The number of satellite 
locations chosen within a given angular sector will be de- 
termined by the spatial resolution of the cloud density and 
intensity required in the model calculation. 

The equations of motion (1) - (3) may be taken to define 
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a set of coupled first order differential equations which 
defines the time evolution of (X/y/Z,x,S’,z) . These equations 
are solved numberically using a fourth order Runge-Kutta 
method. The equations (1) - (3) exhibit no integrals of 
motion. Accuracy of the numerical solution is achieved by 
using carefully constructed sets of small time steps for 
integrating each atom orbit. In perfoming these integrations, 
properly constructed weight factors along each orbit to sim- 
ulate both ionization and D-line excitations of the sodium 
cloud atoms may also be calculated and used to compute the 
column density, D-line intensities, and D-line emission 
profiles of the cloud. 

In the model the ejection of sodium by lo is described 
by atoms emitted radially from a spherical exobase taken to 
lie at an altitude of 780 km above the satellite surface. 

This three dimensional model contains 1298 source points dis- 
tributed uniformly over the exobase surface. Model calcu- 
lations considered in this paper are however restricted to 
monoenergetic ejection of sodium from the inner hemisphere 
of the satellite and involve only 685 of these exobase points. 
The inner hemisphere is selected, since earlier modeling 
efforts (Smyth and McElroy, 1978; Matson et al . , 1978; Macy 
and Trafton, 1980) have shown it to be the dominant emission 
hemisphere to explain the two-dimensional intensity features 
of the sodium cloud. 

Ionizations of the sodium atoms along their orbits in 
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the following model calculations are limited to the simplest 
description of a constant cutoff lifetime of 20 hours and 
no attempt is made, as discussed earlier, to address the more 
complex spatially non-uniform and time-dependent tailoring 
of the sodium cloud produced by its interaction with the lo 
plasma torus. The small lifetime for sodium atoms of order 
one hour in the warm inner region of the lo plasma torus, 
estimated from Voyager plasma data (Bagenal and Sullivan, 
1981) , is in marked contrast to the lifetime estimates of 
15 to 20 hours (Smyth and McElroy, 1978) and 28 hours (Matson 
et al., 1978) obtained from earlier modeling efforts and to 
the 20 hour value to be adopted here. In these earlier 
modeling studies, where no spatial variation of the sodium 
lifetime was assumed, the value of the lifetime was essen- 
tially determined by properly modeling the spatial extent 
of the forward sodivmi cloud observed in Earth-based measure- 
ments. With the presence of the plasma torus since estab- 
lished, it is now clear that the forward extent of the sodium 
cloud is determined by ionization of these cloud atoms upon 
their second encounter with the plasma torus, somewhat ahead 
of lo (as illustrated in Figure 3) , which occurs for flight 
times of order 20 hours. During the first encounter of the 
cloud atoms with the plasma torus (i.e., as they escape lo 
initially) the abundance of sodium atoms in the cloud is 
reduced in a manner which depends critically upon the time 
history of the atom orbits and their location in the plasma 
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torus, which oscillates about the satellite plane. This 
spatial and time dependent destruction of the sodium density 
within the 20 hour cloud envelope will not be treated in this 
paper. The value of 20 hours for the sodium lifetime is 
adopted, however, in order to preserve the correct dimensions 
for the forward sodium cloud. 

From Figure 3, it is relatively easy, however, to anti- 
cipate three effects of including the spatially nonuniform 
ionization of the lo plasma torus on the model calculations 
presented here. Fir.^t note that the forward cloud, as deter- 
mined by the lo plasma torus, will be significantly more elon' 
gated for western phase angles than for eastern phase angles, 
when cloud images are compared for diametrically opposite 
satellite locations about Jupiter. This more elongated be- 
havior of the western cloud has been confirmed in very recent 
measurements of Goldberg et al., (1981) . The actual flight 
times in Figure 3 for which sodium atoms traveling along the 
inner edge of the east and west clouds would have their 
second encounter with the plasma torus are, however, nearly 
equal being about 26 and 24 hours respectively. Second note 
that the density and D-line intensity of the outside portion 
of the sodium cloud will be reduced, relative to the inside 
portion of the sodium cloud, because of its prolonged spatial 
interaction with the lo plasma torus. This may cause the 
east-west cloud distribution asymmetry calculated in this 
paper to be modified in such a way as to reduce the predicted 
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values of the east critical phase angle more than the predicted 
value of the west critical phase angles. The critical phase 
angles are the satellite phase angles for which the sodium 
cloud on the sky plane appears to be distributed approximately 
symmetrically about lo. Third note that because of the first 
ionization encounter of the forward cloud atoms v^ith the lo 
plasma torus, the density and D-line intensities of the sodium 
cloud near lo will be more enhanced, in comparison to the 
forv;ard portion of the sodium cloud, than in the present cal- 
culation. This will cause the east-west intensity asymmetry 
effects that originate near lo to be enhanced and will also 
require that a larger flux of sodixam atoms be emitted from lo 
in order to properly populate the elongated portion of the 
forward cloud. 

Initial velocities of atoms emitted from the staellite 
exobase in the following model calculations are limited to 
values less than 3 km sec , since these lower velocity 
components are those that will contribute dominantly, in the 
actual velocity dispersion, to the cloud intensity for 
brightness levels greater than about 1000 rayleigh (see Smyth 
and McElroy, 19V8y <. It is these larger intensity values v/hich 
are observed in the data of Bergstralh et al. (1975, 1977) 
and Goldberg et al. (1978) . Velocity components greater^ than 
about 3 km sec”^ will, however, contribute significantly to 
the two-dimensional intensity of the sodium cloud for lesser 
brightness levels, such as the peculiar directional features seen 
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by Picher (1980) and Goldberg et al. (1981 ) , and to the 
asynunetric wings of the sodium line profile shape, such as 
has been noted by Trafton (1975 ) , Trauger et al, (1976) and 
Carlson et al. (1978) . 
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3. Model Results: Goldberg's Cloud Distribution Asynvmetry 


The sodium cloud has been observed by Goldberg et al. 
(1978) to exhibit an east-west asymmetry in the spatial dis- 
tribution of the D-line intensities with respect to the phase 
angle location of the satellite. Their observations showed 
that the projected intensity pattern of the cloud on the sky 
plane did not produce a mirror image for satellite orbital 
phase angles separated by 180 degrees. To initially quantify 
and characterize one of the differences observed in compar- 
ison of these east and west cloud images, the critical phase 
angles were selected as parameters. The cloud exhibits two 
critical phase angles which occur for those satellite phc-^se 
angles where the sodium cloud on the sky plane appears nearly 
symmetric about lo as the satellite approaches either the 
eastern or western elongation points of its orbital (i.e., 

90 degrees or 270 degrees, respectively) . 

The first estimates of the eastern and v/estern critical 
phase angles were originally reported by Goldberg et al. 

(1978) as 65 degrees and 230 degrees respectively. These 
values were defined quantitatively by measuring the intensity 
seen in two simulated slits, applied to the two-dimensional 
cloud image and located 9 arc sec on either side of lo, and 
by determining that satellite phase angle for which the ratio 
of the two slit intensities was unity. This method lad to a 
rather well defined value for the critical phase angle in the 
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west, but resulted in a large amount of scatter in determining 
the critical phase angle in the east. The less distinct be- 
havior of the data in the east seemed to reflect a fundamental 
east-west difference in the sodiiim cloud, not a difference in 
the quality of the east and west data. Later estimates made 
with additional observations (Goldberg, 1979) indicated that 
the critical eastern phase angle might even be as large as 70 
to 75 degrees. This then suggested that in the east the central 
axis of the fojrward sodium cloud was inclined only about 15 
to 20 degrees with respect to the line drawn tangent to the 
satellite; orbit at lo, whereas in the west this indlination 
angle was about 40 degrees. The sodium cloud would then have 
an east-west phase lag asymmetry of about 20 to 25 degrees. 

The accuracy of these estimates for the east and west 
critical phase angles and hence for the east-west phase lag 
asymmetry has been, however, somewhat difficult to determine 
with exactness from the data of Goldberg et al. (1978) and 
Goldberg (1979) , since their 2 to 3 hour long measurements 
introduced an orbital smear of 17 to 25 degrees into the 
cloud images. To rectify this situation, Goldberg et al. 

(1981) have recently made several new measurements using an 
improved instrument with an integration time of about 10 
minutes. Their preliminary analysis of these new measurements, 
based upon visual comparisons, indicates that the western 
critical phase angle is distinct and occurs near 235 degrees, 
whereas the eastern critical phase angle is significantly 
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less distinct and occurs somewhere in the angular interval 
from 55 to 70 degrees. The exact value of the east-west phase 
lag asymmetry at present/ therefore/ remains undetermined and 
may or may not be as large as initially suggested. The one 
significant east-west asymmetric cloud feature that does occur 
both in these new measurements and in their older observations 
is the obvious presence of a distinct critical western phase 
angle and the noticeable visual absence of a distinct critical 
eastern phase angle. Quantitative analysis of existing data 
and acquisition of new data will be required to better determine 
the critical eastern phase angle and to understand possible 
variations of its value introduced by the oscillating lo 
plasma torus. Further analysis 6f these data should also 
be useful in identifying and characterizing other salient 
features of the east-west cloud distribution asymmetry/ such 
as the greater elongation of the forward portion of the sod- 
ium cloud observed for western phase angles by Goldberg 
et al. (1981). To provide a step forward/ based upon our 
current understanding, modeling analysis in this section will 
be restricted to exploring and documenting the east-west 
phase lag asymmetry introduced by the solar radiation pressure, 
excluding the spatially non-uniform and time-dependent ioni- 
zation of sodium atoms by the lo plasma torus. 

The preliminary results of calculations in the satellite 
plane, originally present in Paper I and repeated here in 
Figure 1, demonstrate that solar radiation pressure in the ' 
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D-lines causes the envelope of the sodium cloud to behave 
in a manner similar to the observations of Goldberg et al. 

(19 7fi) . Changes in the cloud shape in Figure 1 suggest an 
east -west phase lag of about 15 degrees. It is clear, however, 
that such an estimate is only preliminary and that the actual 
three dimensional distribution of sodium atoms within the 
cloud envelope must be calculated in order to determine a more 
accurate phase lag angle. Suitable three-dimensional model 
calculations for this purpose are presented below. 

Three-dimensional model calculations, illustrating the 
effects of solar radiation acceleration on the D 2 intensity 
distribution of the sodium cloud viewed parallel to the sat- 
ellite plane, are given for satellite phase angles of 55 and 
235 degrees in Figure 4 and Figure 5, respectively. In these 
model calculations, which may be compared directly with the 
simpler results of Figure 1, sodium was emitted radially 

from the inner hemisphere of lo's exobase (2600 km radius) 

—1 

with a velocity of 2.6 km sec and with a cutoff lifetime 

/ 

of 20 hours in the Jovian environment. Each outer contour 

8 “2 

value is 916 Rayleighs for a surface flux of 10 atoms cm 
sec”^, with contours spaced by 367 Rayleighs, increasing 
inward toward lo. The results of Figure 4 and Figure 5 show 
that changes in the D 2 intensity pattern of the sodium cloud 
are not very sensitive to the ill\amination angle of the sun 
relative to the Earth-Jupiter line of sight. The sun angle 
is confined to values between about 167 and 193 degrees as 
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the earth moves about the Sun. For a sun angle of 180 degrees, 
the Earth-Jupiter line of sight and the Sun-Jupiter axis of 
Figure 1 are identical. For a satellite phase angle of 55 
degrees, note that in Figure 4 the effect of solar radiation 
is to rotate the central axis of the forward elongated sodium 
cloud further away from Jupiter as anticipated in the results 
of Figure 1. With solar radiation acceleration included in 
the model calculation, the inner edge of the cloud (right 
portion) in Figure 4 is therefore closer to lo while the outer 
edge (left portion) is extended further from lo. For a phase 
angle of 235 degrees, note that in Figure 5 the solar radiation 
pressure has effectively rotated the sodixam cloud axis more 
toward Jupiter, as anticipated in Figure 1. With solar radi- 
attoi acceleration included, the outer edge of the cloud 
(right portion) in Figure 5 is therefore closer to lo while 
the inner edge of the cloud (left portion) clearly reveals 
that the axis of the forward elongated sodium cloud has al- 
ready begun to swing through our line of sight. 

The model calculations of Figure 4 and Figure 5 indicate 
that solar radiation pressure introduces an east-west phase 
lag asymmetry into the spatial distribution of the sodium 
cloud. This can be illustrated more clearly as follows. In 
the absence of solar radiation acceleration, model calculations 
of the D 2 intensity pattern of the sodium cloud for diametri- 
cally opposite satellite phase angles are presented in Figure 
6, where the sun angle is 180 degrees and the other model 
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parameters and contour levels are the same as in Figure 4 and 
Figure 5. Note that the east and west intensity patterns are 
effectively mirror images of each other, and that the axis of 
the forward elongated cloud can be seen to swing through our 
line of sight at a phase angle of about 55 degrees in the 
east and 235 degrees in the west, yielding no phase lag. 

The observed changes in the intensity pattern of the cloud 
with satellite phase angle are almost completely due to the 
changing geometric viewing perspective of the cloud projected 
onto the plane of the sky. In Figure 7, the D 2 intensity 
patterns of the sodium cloud, where the effects of solar 
radiation acceleration have been included, are presented for 
t'h'^ identical diametrically opposite satellite phase angles, 
model parameters, and contour levels used in Figure 6. Near 
eastern elongation, note that the outer contour of the sodium 
cloud appears symmetric about lo near a satellite phase angle 
of between 60 and 65 degrees, whereas near western elongation 
the cloud appears symmetric at a phase angle of 230 degrees. 
This produces an east-west phase lag asymmetry between 10 and 
15 degrees. The distribution of the D 2 intensity within the 
outer contour in Figure 7 is approximately symmetric for a 
phase angle of 65 degrees in the east and 235 degrees in the 
west. Using this as a criteria, the east-west phase lag asym- 
metry is then about 10 degrees. Note that in Figure 7, the 
passage of the forward cloud axis through our field of view is 
much less distinct in the east than in the west. This mcxkes it 
more difficult to accurately determine the eastern critical 
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satellite phase angle, a persistent characteristic noted in the 

observational data of Goldberg et al. (1978, 1981) . 

The results Figure 6 and Figure 7 are calculated assuming 

a monoenergetic emission velocity of 2.6 km sec""^. Results 
for lower emission velocity components that would be included, 
if a more realistic initial emission velocity distribution 
were adopted, have also been calculated. For emission velo- 
cities of 2.0 km sec”^, 2.2 km sec~^ and 2.4 km sec~^, the 
intensity patterns on the sky plane are again essentially 
mirror images, in the absence of solar radiation pressure, 
when the east and west clouds are compared for diametrically 
opposite satellite phase angles. The critical phase angles 
increase with decreasing emission velocity, having a value, 
in the absence of solar radiation pressure, in the east 
between 55 and 60 degrees for an initial velocity of 2.4 km 

sec”^ and a value of about 60 degrees for initial velocities 

-1 -1 

of both 2.2 km sec and 2.0 km sec . With solar radxation 
pressure included, these model calculations have the same 
general cast-west character as the results presented in Figure 7. 
Based upon the outer contour criteria, the critical eastern 
phase angle is about 65 degrees for the emission velocities 
of 2.4 km sec ^ and below, while the critical western phase 
angle is between 230 and 235 degrees. This provides an 
east-west phase lag asymmetry of about 10-15 degrees. Based 
upon the distribution of intensity within the outer contour, 
the critical eastern phase angle is about 70 degrees while 
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the critical western phase angle is about 235 degrees for 

. -1 
emission velocities of 2.4 km sec ' and lower, producing a 

15 degree east-west phase lag asymmetry. 

-1 

Model results for emission velocities of 2.2 km sec , 

-1 -1 

2.4 km sec and 2.6 km sec are compared in Figure 8 and 
Figure 9 for diametrically opposite phase angles east and 
west of Jupiter. These D 2 intensity images include the 
effects of solar radiation pressure and are calculated for 
the same contour levels, emission conditions, and lifetime 
conditions chosen in Figures 4-7. As the emission velocity 
decreases, the cloud is more closely confined near lo and 
provides a rapidly increasing contribution to the intensity 
of the cloud near the satellite. The exact east and west 
critical phase angles and corresponding east-west phase lag 
asymmetry predicted from these model calculations will then 
depend to some extent upon the shape of the initial velocity 
distribution. Prom the above calculations, however, it is 
clear that a phase lag of about 10 to 15 degrees would be 
present even if an initial velocity distribution were adopted. 
For such an adopted distribution, the passage of the cloud 
through its eastern critical phase angle would also appear 
to occur slower than the passage of the cloud through its 
western critical phase angle. The critical phase angle in 
the east should be within the angular range from 55 to 70 
degrees, while in the west it should be between 230 and 235 
degrees . 
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The presence of .'-i distinct critical western phase angle 
and of a significantly less distinct critical eastern phase 
angle in the observations of Goldberg et al. (1981) suggests 
that the non-uniform ionisation of sodium introduced by the 
lo plasma torus oscillating about the satellite plane may 
tailor and modify, but will not eliminate, the basic charac- 
teristics of the east-west cloud distribution asymmetry de- 
termined by the model calculations in this section. Further 
quantitative analysis of the recently acquired data of Goldberg 
et al. (1981) and future analysis of additional sodium data 
anticipated during next apparition by Goldberg (1981) is 
warranted. These data should provide the necessary information 
for detailed modeling of the east-west cloud distribution 
asymmetry as well as the alternating north-south asymmetry 
discovered by Trafton and Macy (1975) and discussed further 
by Trafton (1980) Improvements in the present sodium cloud 
model to include the non-uniform ionization of the lo plasma 
torus are now in progress. The improved model should provide 
the necessary tool for interpretation of both the east-west 
and the north-south asymmetric phenomena and, thereby, lead 
to a better understanding of the local atmosphere of lo and 
its complex interactions with the Jovian magnetosphere. 
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4. Model Results: Bergstralh's Intensity Asymmetry 

The sodium cloud has been observed by Bergstralh et al. 
(1975,1977) to exhibit an east-west asymmetry in its absolute 
intensity. In their observations, the intensity of the so- 
dium cloud was measured as a function of phase angle through 
a 3" X 8” slit centered on the satellite with the long di- 
mension of the slit oriented approximately perpendicular to 
the projection of lo's orbital plane. This observational 
geometry is illustrated in Figure 10. The long dimension of 
the slit is only slightly larger than the diameter of the 
Lagrange sphere of lo, within which the gravitational field 
of the satellite dominates the planetary gravitational field. 
The spatial region of the cloud sampled by the slit is there- 
fore composed of a mixture of the dense sodium that likely 
exists within the Lagrange sphere and the somewhat less dense 
sodium that occurs, along the line of sight, in the foreground 
and background of this sphere. 

The D-line intensities of the cloud, measured through 
this observing slit, were shown by Bergstralh and colleagues 
to be correlated with lo's orbital position around Jupiter. 
This correlation made it possible for them to identify res- 
onance scattering of sunlight by sodium atoms as the dominant 
excitation mechanism for the cloud. In addition, the measured 
intensity of the cloud with the satellite east of Jupiter was 
shown, in their initial paper (Bergstralh, Matson and Johnson, 
1975) , to be consistently about 50% greater than the measured 
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intensity with the satellite west of Jupiter. In their second 

paper Bergstralh et al,, 1977)/ which included additional ob- 

< 

servations and a presentation of more select data/ this east- 
west intensity asymmetry was more accurately estimated to be 
20 to 25%. 

In paper I, this intensity asymmetry was associated with 

the ability of solar pressure to compress the east cloud and 

expand the west cloud envelope, as illustrated in Figure 1. 

To more fully investigate how solar radiation pressure provides 

a mechanism to explain this intensity asymmetry/ a number of 

three-dimensional model calculations have been performed and 

are presented here. The Dn and D., intensities seen through 

a 3" X 3" slit centered on lo have been calculated as a 

function of both the satellite phase angle and the initial 

speed of sodium atoms emitted from the inner hemisphere of 

Io‘s exobase. In all of these model computations/ the two- 

dimensional distributions of both the column density and the 

intensity of the cloud on the sky plane, as well as the shape 

of the sodium D-line profiles, have also been simultaneously 

calculated, both excluding and including solar radiation 

pressure, to more fully understand the basic phenomena. 

Model calculations for the relative D 2 intensity seen 

through the 3" x 8" viewing slit, taken here to have nominal 

3 4 

dimensions of 9 x 10 km by 2.4 x 10 km, are presented xn 
Figure 11 and compared with the observations of Bergstralh 


28 


et al. (1975) . The relative intensity of the model cal- 
culations is shown as a function of the satellite phase angle 

for four different monenergetic emission velocities spanning 

-1 -1 

the range from 2.0 km sec to 2.6 km sec . In each mono- 

t 

energetic model calculation, the same n\imber of atoms were 
emitted from lo's exobase. These relative intensities were 
calculated for sodium atoms emitted from the inner hemisphere 
of the exobase and include the shadovr effect of the disk of 
lo blocking the cloud behind it. The tilt angle of the 
satellite plane is assumed to be zero relative to the ob- 
server. The sun, observer and Jupiter are assumed to be along 
the same straight line. Intensities calculated for a sa- 

i 

tellite-plane tilt angle of a few degrees, differ only slightly. 

The model calculations of Figure 11 clearly show that 
solar radiation pressure introduces a large east-west inten- 
sity asymmetry for the initial emission velocity of 2.0 km 
sec ^ and a somewhat less pronounced asymmetry for the 2.1 
km sec ^ emission velocity. The asymmetry is present only 
in small amounts for the higher emission velocities of 2.4 

km sec ^ and 2.6 km sec The maximum intensity of the 2.0 

-1 . -1 
km sec results rs 4.38 times larger than the 2.6 km sec 

results in the east and 3.48 times larger in the west. The 

contributions to the intensity for the different velocities 

have a maximum value near 60 degrees in the east and near 

230 degrees in the west. In both the east and west, these 


29 


peak intensities correspond to the enhancement of the sodium 
cloud density produced by the forward portion of the cloud 
as it passes through the line of sight of the viewing slit. 

The peak intensity is larger and more pronounced for the lower 

emission velocities since the forward cloud is then much more 
narrow and is much more tightly confined near lo. This 
spatial behavior of the forward cloud, in the absence of solar 
radiation pressure, was discvsssed earlier by Smyth and McElroy 
(1977) , who at that time noted a similarity between their 
calculated coliamn density peaks and the peaks in the data of 
Bergstralh et al. (1975) . This similarity is also present in 
the results of Figure 11 , although before an accurate compar- 
ison can be made with the observations , the additional modi= 
fications introduced by the spatially non-uniform lifetime of 
the lo plasma torus must be included in the model calculations. 

The east-west ratio of the relative intensity for the 

—1 —1 

2.0 km sec and 2.1 km sec emission velocities is presented 
in Figure 12, where the results of including and excluding 
the effects of lo's shadow are both indicated. For an emission 
velocity of 2.0 km sec ^ the ratio provides a minimum east-west 
asymmetry of 17 percent over most of the angular range when 
the shadow is included, and at a 70 degree phase angle the 
ratio has a peak value of over 1.8. The 2.1 km sec ^ ratio 
is less pronounced, but still provides a significant east-west 
symmetry over much of the phase angle interval. The general 
qualitative behavior of the intensity ratio exhibited in 
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Figure 11 can be simply understood in terms of the actio /,1 of 
solar radiation pressure upon the sodium atom orbits. Thijs 
qualitative explanation will now be discussed in some detail. 

In Figure 12, the intensity ratio is greater than 
unity for east/west phase angles in the range from 0/180 
to 40/220 because the column density of the east cloud in 
the viewing slit is larger than that of the west cloud by 
approximately 20% from an emission velocity of 2.0 km sec . 
The radial velocity v of most sodium cloud atoms with 
respect to the sun is larger, however, in the west cloud 
than in the east cloud, for this phase angle interval, so 
that the D 2 intensity ratio is reduced somewhat from this 
column density ratio through the factor y^(v) in equation 
(14) . This larger velocity of sodiiam atoms in the west cloud 
is also evident in sodium line profiles when model results 
are compared for diametrically opposite satellite phase angles . 
For an east/west phase angle near 50/230 in Figure 12, the 
□2 intensity of the west cloud in the viewing aperature has 
reached its maximum value, because of the phase lag asymmetry 
alignment of the forward cloud. The east cloud, however, 
reaches its maximum D 2 intensity value about 10 degrees later 
and maintains intensity values near this maximum level over 
a longer phase angle interval, beyond the peak, than does 
the west cloud (see the model calculation in Figure 11) . 

This gives rise in Figure 12 to the dramatic increase in the 


intensity ratio within the east/west phase angle range 
from 50/230 to 120/300, and to the asymmetric slope of the 
□2 intensity ratio about its peak value located at a phase 
angle of 70/250. 

For larger east/west phase angles between 120/300 and 
165/345, the D 2 intensity ratio decreases and eventually 
becomes less than unity. In this angular range, most atoms 
now seen in the viewing slit are near lo and not in the 
forward cloud. Solar radiation pressure, in this angular 
range, so acts on these cloud atoms to elevate the column 
density of the west cloud by as much as 18% for the 2.0 
km sec ^ emission velocity and as much as 4% for the 2.1 
km sec ^ emission velocity. In the west cloud, the radial 
velocity v of most atoms with respect to the s\in is also 
enhanced over the east clodd and, through the factor (v) 
in equation (14), causes the D 2 intensity ratio to be re<luced 
ever further. 

The significant enhancement of the calculated column 
density of the east cloud relative to the west cloud within 
the phase angle range 0/180 to 40/220 in Figure 12, may be 
attributed to two factors. The first one is that solar ra- 
diation pressure produces an east cloud that is spatially 
more tightly confined near lo and hence has elevated densities. 
The second and more important factor is that solar radiation 
pressure modifies the escape process of atoms from Xo so that 
the total number of atoms within the east cloud is larger. 
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This modification of the escape process is also largely 
responsible for the elevation of the column densities of the 
west cloud relative to the east cloud within the east/west 
phase angle range form 135/315 to 170/350. 

With regard to the first factor, the density enhancements 

of the east cloud and the density reductions of the west 

cloud within the phase angle range 0/180 to 40/220 in Figure 

12 are more pronounced for atom emission velocities near 
-1 

2.0 km sec . This occurs because many of these atoms, 
shortly after emission from the exobase, experience a near 
balance between the effective gravitational field of Jupiter 
and lo. This near balance causes the relative magnitude of 
the solar radiation acceleration experienced by the sodium 
atom to be favorably amplified. The vector direction of the 
solar radaiation acc»i 1 ,» cation, being oppositely oriented 
in the east and west cloud, then asymmetrically and somev;hat 
abruptly alters the cloud atom orbits. Future evolution of 
these altered atom orbits then produces the density en- 
hancement of the east cloud and density reduction of the 
west cloud. For the higher emission velocities of 2.4 km 
sec"^ and 2.6 km sec”^, this amplification effect is much 
smaller since the sodium atoms, moving more rapidly away 
from lo, spend less time near the Lagrange sphere. The re- 
sulting east-west density changes for these higher emission 
velocities are therefore much less pronounced and result 
primarily from the smooth and continuous action of solar 
radiation pressure over the 20 hour lifetime of the atoms. 
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In the case of the higher emission velocities, the compression 
of the east cloud and the expansion of the west cloud, as 
shown in Figure 1, are, in fact, to first order compensated 
by the column integration along the line of sight and con- 
tribute little to the intensity asymmetry. 

With regard to the second and more important factor, 
the ability of solar radiation pressure to modulate the total 
number of sodium atoms within the cloud occurs, essentially 
for those emission velocities near 2 km sec”^ where the 
escape of atom orbits is energetically possible only from a 
limited surface area of the satellite exobase. This ener- 
getically possible areas of the exobase defines the window 
in the Lagrange sphere surrounding lo through which the atoms 
escape. It is this effective exospheric area, from which 
the emitted atoms are able to escape from the satellite, that 
is modulated by the action of solar radiation. 

The Lagrange window in lo ' s orbital plane for atoms emitted 
from the exobase of the satellite has been discussed, in the 
absence of solar radiation pressure, by Smyth and McElroy (1977) . 
One of their illustrations is shown in Figure 13 for radial 
emission velocities of 1.8 km sec 2.0 km sec and 3.0 km 
sec For a 1.8 km sec ^ velocity, the Lagrange window is 
closed except near the center of the inner satellite hemisphere 
(i.e. at 180®) and near the center of the outer satellite hem- 
isphere (i.e. at 0®). The window is preferentially open in 
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this direction because it is the easiest direction for escape 

of atoms, since the two collinear Lagrange points, and 

nearest to lo lie along the line connecting lo and Jupiter 

(see Figure 2) . For an emission velocity of 2.0 km sec”^, 

the Lagrange window is closed on a small sector of the leading 

inner quadrant and a small sector of the trailing outer quadrant 

-1 

of the exobase. For a velocxty of 3.0 km sec , the Lagrange 
window is completely open, since all atoms emitted radially 
from the exobase now have sufficient energy to escape lo re- 
gardless of their initial surface location. In the absence of 
solar radiation pressure, the percentage of atoms emitted from 
the inner hemisphere of the exobase that do not escape lo is 
determined by model calculations to be 60.3, 28.6, 9.6, 1.8, 

and 0.19 for monoenergetic emissions velocities of 1.8 km sec 

-1 -1 -1 
1.9 km sec , 2.0 km sec , 2.1 km sec-1 and 2.2 km sec 

These percentage values are of course independent of the 

phase angles at which the satellite is observed. 

The ability of solar radiation pressure to modulate, 
as a function of the satellite phase angle, the percentage of 
atoms not escaping from lo, having been initially emitted ra- 
dially from the inner hemisphere of the exobase, is illus- 
trated in Figure 14. Model results, where solar radiation 
pressure has been both included and excluded, are compared 
for monoenergetic emission velocities of both 2.0 km sec*^ 
and 2.1 km sec”*^. The functional behavior exhibited is 
essentially similar in nature for both emission velocities. 
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Solar radiation pressure enhances the escape of sodium atoms 

for smaller satellite phase angles in the east (i.e. between 

0 and 68 degrees for 2,0 km sec”^ and between 0 and 98 degrees 
“1 

for 2.1 km sec ) and for larger phase angles in the west 

rw 1 

(i.e. between 208 and 360 degrees for 2.0 km sec and between 
291 and 360 degrees for 2.1 km sec ^) , This behavior accounts, 
in Figure 12, for the earlier noted enhancement of the column 
density of the east cloud relative to the west cloud in the 
east/west phase angle range of 0/180 to 40/220, and of the 
west cloud relative to the east cloud in the east/west phase 
angle range of 135/315 to 170/350. The general behavior 
exhibited by the effect of solar radiation pressure in Figure 
14 may also be anticipated on geometrical grounds by simply 
considering, as a function of the observed phase angle of the 
satellite, the 20 hour time history of the orientation of 
the solar radiation acceleration vector and that portion of 
the inner hemisphere exobase from which atom escape is mod- 
ulated. 

The maximum observed D 2 intensity reported by Bergstralh 
et al. (1977) using their best data base was 241 k Rayleighs. 
This was measured for a satellite phase angle near eastern 
elongation and assumed that only a 3 x 3 second portion of 
their observing slit was uniformly illuminated. If the 
complete 3x8 arc second slit is considered, this uniform 
intensity is reduced to about 90 k Rayleighs. If this 
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intensity were produced by the maximum D 2 intensity given [ 

by the 2.0 km sec ^ calculation of Figure 11, a satellite 

3 —2 ""X 

flux of order 1 x 10 atoms cm sec would be required, in 
agreement with an earlier estimate by Smyth and McElroy 
(1978). This flxix value of 1 x 10®cm”^sec~^ will, however, 
be revised upward, as noted by Brown (1981) , upon inclusion 
in the model calculation of the highly spatially, non-uniform 
ionization of sodium atoms by the lo plasma torus. 

For the results in Figure 11 and Figure 12 to explain 
the observed east-west intensity asymmetry, the relative 

I 

importance of the emission velocity components near 2.0 km I 

sec”^ must be significant in the actual initial velocity s' 

distribution of sodium atoms in the lo atmosphere. For 

velocity components of 1.8 km sec”^ and lower, ballistic 

orbits are dominant, and the cloud becomes less bright and 

is mostly confined inside the Lagrange sphere, where the 

asymmetric mechanism is not operative. For velocity compo- 

-1 

nents of 2.4 km sec and higher, essentially all the cloud 
atoms escape the Lagrange sphere in a sufficiently short 
time to render the asymmetric mechamism ineffective, and the 
cloud becomes less bright near lo and more spatially extended. 

Inversion of the sodium data of Bergstralh et al. (1975, 1977) 
should therefore be sensitive to the velocity components near 
2.0 km sec”^. This inversion is clearly a desirable modeling 
goal and will be critically dependent upon the calculated results 
presented in Figure 11. The results of Figure 11, however. 
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must first be recalculated so as to include the non-uniform 
spatial lifetime of sodium atoms produced by the presence of 
the lo plasma torus. With these new results, which are 
currently under evaluation, such an inversion process would 
be meaningful. 

Some of the modifications to the results of Figure 11 
and Figure 12 that will be produced by including the presence 
of the lo plasma torus in the new model calculations may be 
anticipated on simple geometrical grounds. The atoms emitted 

from the inner hemisphere of the satellite exobase with a 

-1 

velocity of 2.0 km sec will create a forward cloud that is 
almost completely confined between the inner and outer boun- 
daries of the plasma torus shown in Figure 3. The continual 
exposure of this sodium to the ionization of the plasma torus 
will cause the atom population in the forward cloud to be sig- 
nificantly diminished relative to the atom population near lo. 
As the emission velocity increases, however, less and less 
of the forward cloud will be confined within the plasma torus, 
as illustrated for an emission velocity of 2.6 km sec ^ in 
Figure 3. The ability of the plasma torus to diminish the 
relative number of atoms in the forward cloud will therefore be 
less severe as the emission velocity increases. The relative 
values of the calculated intensity curves in Figure 11 for 
different emission velocities will thus be altered* 

The dramatic increases in the calculated D 2 intensities 
in Figure 11 near 60 degrees and 230 degrees for the emission 
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. . —1 —1 
velocities of 2.0 km sec and 2.1 km sec , which occur 

because of the passage of the forward cloud through the field 

of view of the observing slit, will also be reduced by this 

ionization of the plasma torus. The extent of this reduction 

will depend upon the mean plasma conditions, which have been 

observed to change over a time period of several months (Sandel 

et al., 1979) , and may also depend upon the time modulation 

of these mean conditions by the oscillation of the lo plasma 

torus about the satellite plane. The resulting peak of the 

east/west intensity ratio in Figure 12 at an east/west angle 

of 70/250 will therefore be reduced. The new model predicted 

east-west intensity asymmetry should then be more strongly 

dependent upon the ab.i,lity of solar radiation pressure to 

modulate, as a function of phase angle, the percentage of 

atoms escaping from lo, which was illustrated earlier in 

Figure 14. 
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5 , Concluding Remarks 


A three- dimens ion ptl model for Ig's sodium cloud incor- 
poirtatin'j the full effects of solar radiation acceleration 
associated with solar resonance scattering of atoms in the 
and D 2 lines has been developed. The model has been 
applied to interpret two different east-west asymmetries 
exhibited by the lo sodium cloud in data obtained from 
ground-based telescopes. The major objective of the modeling 
analysis presented was not to invert the sodium data and 
deduce information about the escape mechanism of gases from 
lo. It was rather to provide the preliminary step of docu- 
menting the physical changes produced in the sodivim cloud by 
the solar radiation acceleration, as a function of the satel- 
lite phase angle, and to further clarify, beyond the efforts 
of Paper I, how these changes provide an explanation of the 
observed east-west asymmetries. To provide model results 
appropriate for inversion of the east-west asymmetry data, 
the present calculations would have to be upgraded to include 
the spatially non-uniform and time variable ionization of 
cloud atoms by the lo plasma torus . Such an effort to improve 
the model is in progress. 

Model results for the observed east-west cloud distri- 
bution asymmetry of Goldberg et al, (1978) , Goldberg (1979) 
and Goldberg et al, (1981) were presented in Section 3. The 
observations indicated that the projected intensity patterns 
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of the sodium cloud on the plane of the sky did not produce 
mirror images at orbital phase angl^^S separated by 180 degrees. 
When the satellite was west of Jupiter, the central axis of 
the forward portion of the sodium cloud was observed to pass 
through the line of sight of the observer at a distinct cri- 
tical phase angle between 230 and 235 degrees. When the sat- 
ellite was east of Jupiter, however, the passage of the central 
axis of the forward cloud through the line of sight of the 
observer was significantly slower and less distinct. The 

I 

critical phase angle in the east at present remains uncertain 
within the angular range of 55 to 70 degrees. The east-west 
phase lag asymmetry associated with the east and west critical 
phase angles therefore remains uncertain and may or may not 
be as large as 15 to 20 degrees. Model results, calculated 
for different monenergetic emission velocities, reproduced 
the correct phase angle behavior of the observations, indi- 
cating a distinct critical phase angle in the west between 
230 and 235 degrees and a significantly less distinct critical 
phase angle in the east between 55 and 70 degrees. These 
results, which have excluded the effects of the spatially non- 
uniform ionization of sodium by the plasma torus, suggest an 
east-west phase lag asymmetry of 10 to 15 degrees. 

The behavior of the cloud distribution asymmetry and the 
value of the phase lag asymmetry predicted by the model cal- 
culations are not critically dependent upon the exact distri- 
bution of the emission velocity chosen between 2 km sec ^ 
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and 3 km sec , because these effects depend primarily upon 
the simple geometrical orientation of lo's orbit and the solar 
radiation force vector. In brief/ cloud atoms observed at 
eastern satellite phase angles are initially emitted at 
western phase angles where the solar radiation force acts in 
a direction to enhance either the angular or outward radial 
component of velocity of most atoms over much of their 20 hour 
lifetime. This increases the kinetic energy of these atoms 
relative to Jupiter and produces a cloud that remains more 
distant from Jupiter and closer to lo than would occur 
in the absence of the solar radiation acceleration (see 
Figure 1) . Conversely/ cloud atoms observed at western 
satellite phase angles are initially emitted at eastern phase 
angles where the solar radiation force acts in a direction 
to initially decrease the angular component and later to 
decrease the outward radial component of velocity of most 
atoms during much of their 20 hour lifetime. This initial 
energy loss relative to Jupiter and the later inward directed 
force will produce a cloud which lies closer to Jupiter and 
more distant from lo than would occur in the absence of the 
solar radiation acceleration (see Figure 1) . 

Model results for the observed east-west intensity 
asymmetry of Bergstralh et al. (1975/ 1977) were presented 
in Section 4. The observations indicated that the intensity 
of the cloud seen through a 3 x 8 arc second viewing slit 
centered on lo was about 20 to 25 percent brighter for 
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eastern satellite phase angles than for western satellite 
phase angles. Model calculations exhibited a somewhat sim- 
ilar east~west intensity asymmetry, but only for a narrow 
range of the emission velocities, unlike the model ra^tvilts 
for the east-west cloud distribution asymmetry. Those emission 
velocities which contribute significantly to the intensity 
asymmetry have speeds very near the escape velocity of atoms 
emitted from the satellite exobase. 

The east-west intensity asymmetry predicted by the 
model calculations is critically dependent upon those atoms 
which have emission velocities near the satellite escape speed, 
because these atoms spend a significant amount of time in 
their orbits close to the Lagrange sphere of the satellite, 
where the effective gravitational fields of Jupiter and lo 
are nearly balanced. This near balance amplifies the relative 
importance of the solar radiation force, which then acts in 
two distinct ways to produce the intensity asymmetry. First 
there is a significant enhancement of the normal compression 
of the east cloud and the expansion of the west cloud, illus- 
trated in Figure 1, and an additional geometrical amplification 
provided by the phase lag asymmetry. Second and more impor- 
tantly, for a constant emission flux from the exobase, there 
is a periodic modulation with respect to the satellite phase 
angle of the total number of sodium atoms that escape the 
satellite, and hence, that populate the extended sodium cloud. 
The combination of these enhancement and modulation processes 
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produces an elevated D 2 intensity in the east cloud relative 
to the west cloud over a significant east-west phase angle 
interval . 

In sununary, the model predicted changes in the shape and 
spatial orientation of the sodium cloud envelope and in the 
intensity distributions within the envelope, provide a simple 
explanation for the east-west oloud distribution asymmetry 
and for the east-west intensity asymmetry. These results 
follow directly from the lack of circular symmetry about Ju- 
piter of the forces acting on the cloud atoms. The lack of 
this circular symmetry is introduced by the presence of the 
solar radiation force in the model calculations. 

The identification of the east-west asymmetry in the 
model calculations with sodium atoms emitted within a some- 
what narrow range of the velocities near the satelllite escape 
speed i.s significant. It both provides an explicit and simple 
mechanism to explain the observed east-west intensity asymmetry, 
and suggests that if there are no other major contributors to 
the intensity asymmetry, this narrow portion of the velocity 
dispersion must be well represented in the actual initial 
velocity distribution of sodium atoms operative in the lo 
atmosphere. Careful model inversion of the east-west inten- 
sity data would then be sensitive to accurately determining 
the relative weights for the contributing velocity components 
in the 1.8 km sec ^ to 2.2 km sec ^ velocity interval. This 
analysis would be somewhat complementary to the inversion of 
the two-dimensional intensity data for the sodium cloud on the 
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sky plane (such as the phase lag asymmetry data) , which is 
primarily sensitive to the emission velocities between 2 km 
sec ^ and 3 km sec and to the line profile data of the 
sodiuKt cloud, which is particularly sensitive to higher 
emission velocities, up to 15 km sec”^ or 20 km sec”^. A 
consistent and simultaneous analysis of all three types of 
sodixam data should provide the necessary overlapping constraints 
to deduce the real initial velocity distribution function for 
sodium atoms escaping from lo. 
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FIGURE CAPTIONS 


Figure 1 


Figure 2 


Figure 3 


East-West Changes in the Sodium Cloud. Calculations 
performed in Paper I for the solar radiation per- 
turbed cloud shape, indicated by the shaded area, 
and for the unperturbed cloud shape, indicated by 
the dashed boundary, are compared for diametrically 
opposite satellite phase angles. 

Dimensionless Coordinate System for the Modified 
Circular Restricted Three Body Problem in lo's 
Orbital Plane. The locations of the five classical 
Lagrange points at L^ , and , the 

satellite at x^, and the planet at in the (x,y)- 
plane are indicated. 

Interaction of the Sodium Cloud and the Plasma 
Torus. The spatially projected overlap in the 
satellite plane of the lo plasma torus with the 
solar radiation perturbed cloud shape (solid boun- 
dary with heavy shading) and with the unperturbed 
cloud shape (heavy dashed boundary) of Figure 1, 
is indicated for the diametrically opposite sa- 
tellite phase angles. 


f 


48 



Figure 4 


Effects of Solar Radiation Pressure in the East 


k 


Cloud. Calculated D 2 intensity contours of the 
sodium cloud are compared, with and without the 
effects of solar radiation pressure, for the in- 
dicated phase and sun angles. The outer dashed 
lines are for reference in comparison of the left 
and right cloud boundaries, and th-fs cer^ter dashed 
line is for vertical alignment of the satellite 
location. 

Figure 5. Effects of Solar Radiation Pressure in the West 
Cloud. See legend to Fig. 4. 

Figure 6. East-West Comparison of the Sodium Cloud. Model 
calculations for the D 2 intensity contours of the 
sodiiam cloud, excluding the effects of solar ra- 
diation pressure, are shown for diametrically 
opposite satellite phase angles east and west of 
Jupiter. No east-west phase lag asymmetry is 
evident. See text for discussion. 

Figure 7. East-West Comparison of the Sodium Cloud. Model 
calculations for the intensity contours of the 
sodium cloud, including the effects of solar ra- 
diation pressure, are shown for diametrically op- 
posite satellite phase angles east and west of 
Jupiter. A distinct east-west phase lag asymmetry 
is evident. See text for discussion. 
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Dependence of the Sodium Cloud on Emission Velocity. 


Figure 8, 

Model calculated D 2 intensity contours for the 
sodium cloud, including solar radiation pressure 
effects, are shown for diametrically opposite sa- 
tellite phase angles and compared for emission 

. -1 -1 
velocxtxes of 2.2 km sec , 2.4 km sec and 2,6 

km sec”^. See text for discussion. 

Figure 9. Dependence of the Sodium Cloud on Emission Velocity. 
See legend to Pig. 8. 

Figure 10. Observational Slit for the Intensity Asymmetry 

Measurements. The rectangular observational slit 
of Bergstralh et al. (1975, 1977) , centered on lo, 
is shown in relation to the sodium cloud, the 
satellite orbit and the planet. 

Figure 11. Comparison of Observation and Model Results for 

the East-West Intensity Asymmetry. The observations 
of Bergstralh et al. (1975) are compared to model 
calculations for the relative D 2 intensity seen 
through their rectangular viewing slit. Model 
calculations for four different monoenergetic 
values of the emission velocity are shown. See 
text for discussion. 
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Figure 12. Model Results for the , ' •vSt’^West Intensity Ratio. 

The east to west D 2 intensity ratio of Figure 11 
is shown as a function of the east to west phase 
angles for emission velocities of 2.0 km sec”^ 
and 2.1 km sec”^. 

Figure 13. Lagrange Window in lo's Orbital Plane for Atoms 
Emitted from the Exobase of the Satellite. lo, 
surrounded by its exobase (the dashed circle ) , is 
shown as depicted in an original illustration of 
Smyth and McElroy (1977) . The radial coordinate 
indicates the flight time of an atom emitted ra- 
dially from the exosphere at the angular orien- 
tation about the satellite shown. Upon capture 
by lo, the 5® angular-resolution radial sector 
containing that orbit is blackened for all later 
time. Ballistic orbits are indicated for a capture 
time of 50 hr. For greater capture times, black- 
ening indicates that atoms, initially on an escape 
trajectory, have been recaptured because of a close 
encounter with lo. 
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Figure 14 . Solar Radiation Pressure Modulation of the Sodium 
Atom Escaping from lo. The percentage of sodiimi 
atoms not escaping lo, which are emitted uniformly 
and at a constant rate from the inner hemisphere 
of the satellite exobase with monoenergetic 
emission velocities of 2.0 km sec ^ and 2.1 km 
sec are compared as a function of the lo phase 
angle for model calculation excluding (dashed 
curve) and including (solid curve) the effects 
of solar radiation pressure. See text for dis- 
cussion. 
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Figure 11 
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Figure 13 



